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Three  computer  programs  were  written  with  the  objective  of  predicting 
the  structural  weight  of  aircraft  through  analytical  methods.  The  first 
program,  the  structural  weight  estimation  program  (SWEEP) , is  a completely 
integrated  program  including  routines  for  airloads,  loads  spectra,  skin  tem- 
peratures, material  properties,  flutter  stiffness  requirements,  fatigue  life, 
structural  sizing*  and  for  weight  estimation  of  each  of  the  major  aircraft 
structural  continents.  The  program  produces  first-order  weight  estimates  / / 
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and  indicates  trends  when  parameters  are  varied.  Fighters,  bombers , and  cargo 
aircraft  can  be  analyzed  by  the  program.  The  program  operates  within  100,000 
octal  units  on  the  Control  Data  Corporation  6600  computer.  Two  stand-alone  pro- 
grams operating  within  100,000  octal  units  were  also  developed  to  provide 
optional  data  sources  for  SWEEP.  These  include  (1)  the  flexible  airloads 
program  to  assess  the  effects  of  flexibility  on  lifting  surface  airloads,  and 
(2)  the  flutter  optimization  program  to  optimize  the  stiffness  distribution 
required  for  lifting  surface  flutter  prevention. 

The  final  report  is  composed  of  11  volunes.  This  volume  (Volume  XI)  con- 
tains the  methodology,  program  description,  and  user's  information  for  the 
flexible  loads  stand-alone  program. 
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“R 


Sv 
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CLA 


CYB 

DF  or  6p 


Wing  flap  chord,  in. 

Lifting  surface  lift  coefficient 

Section  lift  coefficient 

Flexible  wing  lift-curve  slope,  per  radian 
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DXW(B)N 


X-distance  from  wing  apex  to  center  of  pressure  of 
unit  airload  on  exposed  wing  for  distribution  N,  in. 


DXWN 
DXW1 
DXW2 
DXW3 

DXW4 
El 

FCJ.D 
F(J,2) 
F(J,3) 
F(J,4) 
FS 

fsea-o 

! 

I 
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X-distance  from  wing  apex  to  center  of  pressure  of 
unit  airload  on  wing,  including  carry-over  load  for 
distribution  N,  in. 

X-distance  from  wing  apex  to  center  of  pressure  of 
airload  on  wing,  including  carry-over  load  for 
flexible  wing  lift  due  to  angle  of  attack,  in. 
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X-distance  from  wing  apex  to  center  of  pressure  of 
airload  on  wing,  including  carry-over  load  for 
flexible  wing  lift  due  to  vertical  acceleration,  in. 

2 

Wing  bending  stiffness,  lb-in. 

Sumnation  of  rigid  wing  unit  loads  on  strips  1 
through  10  for  distribution  J 

Summation  of  flexible  wing  unit  loads  on  strips  1 
through  10  for  distribution  J 

Summation.  of  incremental  flexible  wing  unit  airloads 
on  strips  1 through  10  for  distribution  J 

Ratio  of  flexible  load  to  rigid  load, 

F(J,2)  - (F(J,1) 

Fuselage  station,  in. 

Fuselage  station  at  intersection  of  wing  elastic 
with  plane  of  symmetry,  in. 
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Wing  torsional  stiffness,  lb-in. 

Integer  used  to  designate  wing  strip  number 

2 

Airplanw  pitching  moment  of  inertia,  slug  ft 


Airplane  yawing  moment  of  inertia,  slug  ft 

Integer  used  to  designated  wing -loading 
distribution,  i.e., 

J-l,  loading  due  to  angle  of  attack 
J«2,  loading  due  to  flap  deflection 
J-3,  loading  due  to  vertical  acceleration 

Flap  span  normalizing  parameter 

Flap  lift  effectiveness  parameter, 

(dcj/d 6 p)  t (dCj/da) 

Gust  alleviation  factor 


Span-loading  normalizing  parameter 

Body  nose  length  (x -distance  from  body  nose  to 
cross  section  of  maximum  area),  in. 

Mach  number 

Cotangent  of  leading  edge  sweep  angle 

Bending  moment  along  load  reference  line  in  swept 
reference  system,  in. -lb 

Torsional  moment  along  load  reference  line  about 
load  reference  line,  in. -lb 

Exposed  lifting  surface  panel  rolling  moment  at  body 
interface  station  (for  wing,  horizontal  tail,  or 
vertical  tail),  in. -lb 
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My 

y(SOB) 

Exposed  lifting  surface  panel  pitching  moment  at 
intersection  of  load  reference  line  and  body  inter- 
face biation  (for  wing  or  horizontal  tail),  in. -lb 

Mz 

(SOB) 

Exposed  lifting  surface  panel  yawing  moment  at 
intersection  of  load  reference  line  and  body 
interface  station  (for  vertical  tail),  in. -lb 

N 

Integer  used  to  designate  type  of  wing  airload 
distribution,  i.e., 

N-l,  flexible  lift  due  to  angle  of  attack 
N»2,  rigid  lift  due  to  flap  deflection 
N=3,  incremental  flexible  lift  due  to  flap 
deflection 

N"4,  flexible  lift  due  to  vertical  acceleration 

Ny  or  NY 

Airplane  side  load  factor 

Nz  or  NZ 

Airplane  vertical  acceleration 

nzl 

Airplane  design  limit  vertical  load  factor 

nzm 

Airplane  maneuver  vertical  load  factor 

• 

p 

2 

Airplane  pitching  acceleration,  rad/sec 

P A FLEX 

Incremental  airload  due  to  flexibility 

p 

FLEX 

Flexible  airload  (rigid  load  plus  incremental  load 
due  to  flexibility) 

PRIGID 

Rigid  load 

pwo 

Initial  exposed  wing  panel  load 

PW1 

Initial  wing  panel  load,  including  carry-over  load 
on  body 
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Pyvr  or  PYV  Vertical  tail  airload,  including  carry-over  load,  lb 

Pw  or  PYN  Side  load  on  body  nose,  lb 


PZH  or  PZHB 


PZN  0R  PZN 


Total  vertical  load  on  horizontal  tail,  including 
carry-over  load  on  body,  lb 

Vertical  load  on  body  nose,  lb 


PZW1  or  PZWB1 
PZW2  or  PZWB2 
PZW3  or  PZWB3 
PZW4  or  PZWB4 


Total  wing  airload,  including  carry-over  load  on 
body  for  flexible  lift  due  to  angle  of  attack,  lb 

Total  rigid  wing  airload  due  to  wing  flap 
deflection,  lb 

Total  incremental  flexible  wing  airload  due  to  wing 
flap  deflection,  lb 

Total  flexibility  wing  airload  due  to  vertical 
acceleration,  lb 


AP(ll)  Wing  unit  carry-over  load  on  body,  per  side 

APRoCI)  Unit  airload  on  wing  strip,  I,  for  rigid  wing 

airload  due  to  angle  of  attack 


Unit  airload  on  wing  strip,  I,  for  rigid  wing 
airload  due  to  wing  flap  deflection 


Unit  inertia  load  on  wing  strip,  I,  normalized  to 
an  exposed  wing  panel  weight  of  -1.0  lb 


APZH 

APZW1 
Q or  QDOT 

q 


Increment  horizontal  tail  airload  due  to  gust  or 
pitching  acceleration,  lb 

Increment  wing  airload  due  to  gust,  lb 

2 

Airplane  pitching  acceleration,  rad/sec 
Free -stream  dynamic  pressure,  psf 
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R or  ROOT 


Airplane  yawing  acceleration,  rad/sec 


2 


RLOADS 

"n 


Rigid  loads  at  structural  influence  coefficient 
points 

Body  nose  effective  radius,  in.  ( » >/  (A^j^y/2  n ) 


or  ^ 


Area  of  theoretical  horizontal  tail,  ft 


2 


SlTorSV 


Area  of  theoretical  vertical  tail,  ft 


2 


Area  of  the  theoretical  wing,  ft 


2 


SIC(x,  y ; f,  n) 


[SIC] 

[SIC] 


Structural  influence  coefficient;  i.e.,  deflection 
at  point  x,y  due  to  a 1 -pound  load  at  point  J,  rj-, 
in. /lb 

Matrix  of  structural  influence  coefficients 

Matrix  of  streamwise  slopes 

Shear  along  load  reference  line,  lb 


Sz  Exposed  lifting  surface  shear  at  intersection  of 

(SOB)  load  reference  line  and  body  interface  station,  lb 

UL(I)  Average  running  unit  normal  loading  for  wing 

strip  I 

ULn  Running  unit  normal  load  at  n -station  along  load 

reference  line 


UT(I)  Average  running  unit  torsional  moment  about  load 

reference  line  for  strip  I 

UTtl  Running  unit  torsional  moment  about  load  reference 

line  at  n -station  along  load  reference  line 
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ULN 

UTN 

USZN 

UMXTN 

UMYTN 


UPZWN 

UMYWN 

USZW(B)N 

UMXW(B)N 

UMYW(B)N 

UPZB(W)N 

l?WB(W)N 


Same  as  ULn  or  UL(I),  but  is  normalized  for  a 
1-pound  wing  panel  load,  including  body  carry-over 
load  for  distribution  N,  lb-in. 

Same  as  UT^  or  UTCI)  • but  is  normalized  for  a 
1 -pound  wing  panel  load,  including  body  carry-over 
load  for  distribution  N,  in. -lb-in. 

Unit  normalized  shear  at  stations  along  load 
reference  line  for  distribution  N,  lb 

Unit  normalized  bending  moment  at  stations  along 
load  reference  line  for  distribution  N,  in. -lb 

Unit  normalized  torsional  moment  about  load  reference 
line  at  stations  along  load  reference  line  for 
distribution  N,  in. -lb 

("1.0)  Unit  normalized  wing  panel  load,  including 
carry-over  load  on  body  for  distribution  N,  lb 

Unit  normalized  wing  panel  total  pitching  moment 
(including  carry-over  load)  at  intersection  of  load 
reference  line  and  plane  of  symmetry  for 
distribution  N,  in. -lb 

Exposed  wing  panel  unit  normalized  load  for 
distribution  N,  lb 

Exposed  wing  panel  unit  normalized  rolling  moment 
at  body  interface  station  for  distribution  N,  in. -lb 

Exposed  wing  panel  unit  normalized  pitching  moment 
at  intersection  of  load  reference  line  and  body 
interface  station  for  distribution  N,  in. -lb 

Unit  normalized  carry-over  load  on  body  per  side 
for  distribution  N,  lb 

Unit  normalized  pitching  moment  about  intersection 
of  load  reference  line  and  plane  of  symmetry  for 
body  carry-over  load  per  side  for  distribution  N, 
in. -lb 
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USZW2 


SOB 


Unit  normalized  exposed  wing  rigid  load  due  to  flap 
deflection  («  USZW(B)N  for  N-2),  lb 


Equivalent  airspeed,  knots 


'N 


Body  nose  volume  (over  length,  £N),  in. 


W Airplane  gross  weight,  lb 

Wg(I)  Weight  of  wing  strip,  I,  lb 


Exposed  wing  panel  weight,  lb 


X 

x 

(X/C) 

(X/Oo 

(x/c)F 

XBW(B) 

XBB(W) 

XBH 

XBV 


Distance  measured  along  a line  paralled  to  airplane 
X-axis  or  fuselage  reference  line,  in. 

Chordwise  location  of  wing  structure  influence 
coefficient  point  where  deflection  is  calculated 
(see  Figure  3),  in. 

Center  of  pressure  location  in  fraction  of  chord 

Center  of  pressure  location  in  fraction  of  local 
chord  for  lift  due  to  angle  of  attack 

Center  of  pressure  location  in  fraction  of  local 
chord  for  rigid  lift  due  to  wing  flap  deflection 

Fuselage  station  of  center  of  pressure  of  exposed 
wing  airload,  in. 

Fuselage  station  of  center  of  pressure  of  wing 
carry-over  airload  on  body,  in. 

Fuselage  station  of  center  of  pressure  of 
horizontal  tail  airload,  including  carry-over 
load,  in. 

Fuselage  station  of  center  of  pressure  of  vertical 
tail  airload,  including  carry-over  load,  in. 
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AXA  Distance  aft  of  load  reference  line  measured  normal 

to  load  reference  line  (see  Figure  3),  in. 

(X/C)  _ Location  of  aft  influence  coefficient  points  in 

fraction  of  chord  (see  Figure  3) 

^C^FWD  Location  of  forward  influence  coefficient  points  in 

fraction  of  chord  (see  Figure  3) 


(X/C1  Location  of  center  of  gravity  at  wing  strip  weight 

in  fraction  of  wing  strip  mean  chord  (see  Figure  2) 

Fuselage  station  of  loading  edge  of  wing  local 
chord,  in. 


XCG 

Xjj  or  XBN 


Fuselage  station  of  airplane  center  of  gravity,  in. 

Fuselage  station  of  center  of  pressure  of  body 
nose  load,  in. 

Fuselage  station  of  body  nose  station,  in. 


*RW  or  *WE 


Fuselage  station  of  theoretical  wing  root  chord 
leading  edge,  in. 


Xj^  or  Fuselage  station  of  theoretical  horizontal  tail 

root  chord  leading  edge,  in. 


*RV  or 


Fuselage  station  of  theoretical  vertical  tail  root 
chord  leading  edge,  in. 


XW1 


XW2 


XW3 


Fuselage  station  of  center  of  pressure  of  flexible 
wing  airload,  including  carry-over  load,  due  to 
angle  of  attack,  in. 

Fuselage  station  of  center  of  pressure  of  rigid 
wing  airload,  including  carry-over  load,  due  to 
flap  deflection,  in. 

Fuselage  station  of  center  of  pressure  of  incre- 
mental flexible  wing  airload,  including  carry-over 
load,  due  to  flap  deflection,  in. 
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XW4  Fuselage  station  of  center  of  pressure  of  flexible 

wing  airload,  including  carry-over  load,  due  to 
vertical  acceleration,  in. 

y Spanwise  location  of  wing  structural  influence 

coefficient  point  where  deflection  is  calculated 
(see  Figure  3),  measured  along  elastic  axis,  in. 

Ypj  Spanwise  station  (butt  line)  of  inboard  end  of 

wing  flap,  in. 

YpQ  Spanwise  station  (butt  line)  of  outboard  end  of 

wing  flap,  in. 

Yfij,  Ygg,  YfiW  or  Yj^  Spanwise  station  (butt  line)  of  wing -body  or 

horizontal  tail-body  interface  station,  in. 

Spanwise  station  (butt  line)  of  center  of  pressure 
of  exposed  wing  airload  due  to  distribution  N,  in. 

Spanwise  station  along  load  reference  line  measured 
from  intersection  of  load  reference  line  and  plane 
of  symmetry,  in. 

Spanwise  station  (butt  line)  of  center  of  pressure 
of  exposed  wing  airload,  in. 

Spanwise  station  (butt  line)  of  center  of  pressure 
of  horizontal  airload,  including  carry-over  load, 
in. 

Z-distance  from  vertical  tail  theoretical  root 
chord  to  vertical  tail -body  interface  station,  in. 

Distance  along  vertical  tail  load  reference  line 
measure  from  intersection  of  load  reference  line 
and  theoretical  vertical  tail  root  chord  plane,  in. 

Z-distance  from  theoretical  vertical  tail  root 
chord  to  center  of  pressure  of  vertical  tail 
airload,  in. 


YW(B)N 

Ya 

YBW(B) 

YBH  or  YH 

hi 

Za 

ZBV 
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Angle  of  attack,  radians 
Load  due  to  angle  of  attack 
Vortex  strength,  (C  cp  V/2 

Nondimens ional  span  station,  Y(b/2) 

Nondimens ional  span  station  at  center  of  wing 
strip,  I 

Distance  along  wing  elastic  axis  from  plane  of 
symmetry  to  influence  coefficient  point  at  which 
unit  load  is  applied  (see  Figure  3),  in. 

Sweep  angle  of  wing  elastic  axis  (see  Figure  3) , deg 
Sweep  angle  of  leading  edge,  deg 


Sweep  angle  of  load  reference  line,  deg 

Taper  ratio  of  theoretical  lifting  surface 

Spanwise  distance  along  wing  elastic  axis  measured 
from  plane  of  symmetry  (see  Figure  3) , in. 

Airplane  mass  ratio,  2(W/S)  -r  (Pn  CAV  cLQrWp) 

Air  density,  slug /ft3 

Chordwise  location  of  structural  influence 
coefficient  at  which  unit  load  is  applied  (see 
Figure  3) , in. 


Section  I 


INTRODUCTION 


During  accelerated  flight  conditions  at  high  speeds,  the  deflections  of 
the  lifting  surface  structure  tend  to  redistribute  the  airloads.  The  result- 
ing airload  distribution  can  be  considerably  different  from  that  computed  on 
the  assumption  of  complete  rigidity.  The  redistribution  results  from  the 
change  in  strearawise  angle -of- attack  along  the  span  caused  by  the  torsional 
and  bending  deflections.  For  a given  mach  number,  the  greater  the  dynamic 
pressure,  the  greater  will  be  the  load  redistribution. 

The  airloads  module  in  the  .SWEEP  program  does  not  include  the  effects  of 
aeroelasticity;  i.e. , the  changes  in  airload  distributions  due  to  structural 
deflections.  A significant  refinement  is  obtained  by  including  the  added 
effects  on  loads  caused  by  wing  structural  flexibility.  This  is  accomplished 
by  the  use  of  the  stand-alone  flexible  airloads  program  described  herein. 
Methods  and  formulation  employed  in  the  stand-alone  program  are  presented  in 
Section  II.  The  computer  program  description  and  program  usage  information 
are  presented  in  Sections  III  and  IV,  respectively. 

The  stand-alone  program  requires  a substantial  amount  of  external  input 
data.  These  data  consist  of  (1)  airplane  geometry  data  identical  to  that 
used  by  the  airloads  module  (BLCNTL)  in  the  SWEEP  program,  (2)  the  wing  El 
and  GJ  distribution  and  elastic  axis  location,  and  (3)  the  specific  flight 
condition  case  data.  The  specific  flight  condition  case  data  include  type 
of  flight  condition  (balanced  maneuver,  vertical  or  lateral  gust,  and  pitch- 
ing or  yawing  acceleration) , mach  number  and  altitude  combinations,  limit 
maneuver  load  factors,  pitching  and  yawing  accelerations,  airplane  weight  and 
CG  location,  and  estimated  wing  weight  distribution.  The  program  calculates 
the  airload  and  center- of -pres sure  location  for  each  airplane  component  and 
the  airload  shear,  bending  moment,  and  torsion  distribution  on  the  wing  and 
empennage  surfaces,  all  for  the  specified  flight  condition  case. 
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Section  II 


METHODS  AND  FORMULATION 
FLEXIBLE  AIRLOADS  PROGRAM  FUNCTIONS 


The  objective  of  the  flexible  airloads  stand-alone  program,  BFCNTL,  is 
to  deteimine  the  airloads  on  the  airplane  components,  including  the  effects 
of  wing  flexibility.  These  loads  are  determined  for  a specific  flight  condi- 
tion case  and  are  used  as  an  optional  external  input  to  the  SWEEP  program. 

The  methods  employed  are  described  in  the  order  that  the  subroutines  USPANF, 
BNLDSF,  and  SPAEMF  are  used  in  the  stand-alone  program. 

Subroutine  USPANF  is  used  to  determine  the  lifting  surface  unit  airload 
shears,  moments,  centers  of  pressure,  and  lift  curve  slopes,  all  for  a 
specified  mach  lumber  and  altitude.  For  the  empennage  (horizontal  and 
vertical  tail)  surfaces,  the  methods  described  in  Volune  III,  Section  II  are 
used.  However,  for  the  wing,  the  methods  are  revised  to  include  the  effects 
of  wing  flexibility. 

Subroutine  BNLDSF  is  used  to  determine  gross  limit  airload  and  center  of 
pressure  on  each  of  the  airplane  components  and  the  airplane  inertia  factors 
for  specified  flight  conditions.  The  methods  are  revised  to  include  effects 
of  wing  flexibility. 

Subroutine  SPAEMF  is  used  to  determine  limit  airload  shear,  bending 
moment,  and  torsion  distributions  on  the  lifting  surfaces  for  a specific 
flight  condition.  For  the  empennage  surfaces,  the  methods  described  in 
Volume  III,  Section  II,  are  used.  However,  for  the  wing,  the  methods  are 
revised  to  include  the  effects  of  wing  flexibility. 


BASIC  FLIGHT  CONDITIONS 


The  flexible  airload  stand-alone  program  is  designed  to  calculate 
airplane  structural  component  airloads  for  a specific  flight  condition  case. 
The  specific  flight  condition  case  is  defined  by  the  types  of  condition 
(maneuver  or  gust),  mach  No.,  altitude,  gross  weight,  eg  position,  and  wing 
weight  distribution.  Cases  from  which  a specific  case  can  be  selected  are 
shown  in  Figure  1.  It  is  noted  that  a specific  case  may  consist  of  more 
than  one  type  of  condition  and  at  more  than  one  altitude.  Cases  available 
are  similar  to  those  in  the  SWEEP  program  (Volume  III,  Section  II),  except 
that  additional  altitudes  are  provided  because  effects  of  flexibility  vary 
with  dynamic  pressure  at  a given  mach  lumber. 
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balanced  maneuver 


»pnJHlV 
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Figure  1.  Basic  flight  condition 


Condition  case  No.  1 through  7 and  10  and  11  cover  the  basic  flight 
conditions  with  wing  flaps  retracted.  Condition  case  No.  8 is  a limit  load 
factor  maneuver  condition  with  wing  flaps  down.  Condition  case  No.  9 is  a 
1.0  g flaps-down  trinmed  condition  for  use  with  landing  conditions. 


DETERMINATION  OF  ' IFTING  SURFACE  UNIT  AIRLOADS 

Methods  of  analysis  used  to  develop  lifting  surface  unit  airload  distri- 
bution, surface  lift  curve  slopes,  and  surface  airload  centers  of  pressure 
are  presented  in  this  section.  Unit  airloads  are  defined,  and  basic  data 
used  for  the  determination  of  the  surface  unit  airloads  are  also  presented. 

Unit  airload  shears  and  moments  are  determined  at  13  selected  spanwise 
stations  along  the  selected  load  reference  line,  including  root  and  tip 
stations.  Unit  airload  shear  and  moments  at  the  surface-body  interface 
station  are  determined  in  the  unswept  (body  axes)  system.  Overall  centers 
of  pressure  locations  of  the  exposed  panel  and  body  carry-over  loads  are 
determined  with  respect  to  the  theoretical  surface  apex. 


DEFINITION  OF  UNIT  AIRLOADS 

The  wing  unit  airload  distributions  are  defined  in  the  following 
paragraphs. 


Unit  Rigid  Lift  Due  to  Angle  of  Attack 

Unit  rigid  lift  due  t'>  angle  of  attack  is  the  airload  distribution  due  to 
angle  of  attack  for  a rigid  lifting  surface  load  of  1 pound  per  side  or  panel 
at  a specific  mach  number. 


Unit  Flexible  Wing  Lift  Due  to  Angle  of  Attack 


The  unit  flexible  wing  lift  due  to  angle  of  attack  is  the  net  airload 
distribution  due  to  angle  of  attack  for  the  flexible  wing  surface,  and  is 
normalized  to  a wing  load  of  1.0  pound  per  side  at  a specific  combination  of 
mach  number  and  dynamic  pressure.  The  net  airload  distribution  includes  the 
combination  of  the  rigid  and  the  aeroelastic  increment. 
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Unit  Rigid  Lift  Due  to  Wing  Flap  Deflection 


Unit  rigid  lift  due  to  wing  flap  deflection  is  the  airload  distribution 
due  to  flap  deflection  for  a rigid  wing  surface  load  of  1 pound  per  side. 


Unit  Incremental  Flexible  Lift  Due  to  Wing  Flap  Deflection 

Unit  incremental  flexible  lift  due  to  wing  flap  deflection  is  the 
incremental  flexible  airload  distribution  due  to  the  application  of  the  rigid 
airload  distribution  due  to  wing  flap  deflection,  and  is  normalized  to 
1 pound  per  side  at  a specific  combination  of  mach  number  and  dynamic 
pressure. 


Unit  Rigid  Load  Due  to  Vertical  Acceleration 


Unit  rigid  load  due  to  vertical  acceleration  is  the  wing  weight  dis- 
tribution normalized  to  an  exposed  wingload  of  -1  pound  per  side. 


Unit  Flexible  Wing  Lift  Due  to  Vertical  Acceleration 


Unit  flexible  wing  lift  due  to  vertical  acceleration  is  the  incre- 
mental flexible  airload  distribution  due  to  the  application  of  the  wing 
weight  distribution,  and  is  normalized  to  1 pound  per  side  at  a specific 
combination  of  mach  number  and  dynamic  pressure. 


BASIC  DATA  FOR  UNIT  AIRLOAD  DETERMINATION 

The  basic  data  required  for  unit  ai. loads  include  all  of  the  lifting 
surface  aerodynamic  and  geometry  data  described  in  Volume  III,  Section  II, 
and  the  following  wing  data: 

1.  Exposed  wing  panel  weight  distribution  consisting  of  weight  and 
CG  for  each  of  10  equally  spaced  chordwise  strips,  as  shown  in 
Figure  2.  CG  is  in  terms  of  X/C  of  the  strip  mean  chord. 

2.  Sweep  angle  of  the  wing  elastic  axis,  and  the  fuselage  station 
of  the  wing  elastic  axis  at  the  centerline  of  synmetry, 

FSfiA  = 0,  are  shown  in  Figure  3. 
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3.  The  variation  of  the  exposed  wing  stiffness  parameters,  El  and  GI, 
with  span  station,  yA  , along  the  elastic  axis. 

WING  UNIT  AIRLOAD  DISTRIBWIONS 

The  general  procedure  for  the  calculation  of  unit  flexible  wingload 
data  consists  of  the  following  steps: 

1.  Determination  of  unit  rigid  loads  and  centers  of  pressure  (CP's)  on 
10  equal-span  chordwise  strips  between  the  side  of  the  body  and  wing- 
tip,  and  on  the  strip  between  the  side  of  the  body  and  airplane 
centerline  of  symmetry. 

2.  Determination  of  win#  rigid  lift-curve  slope. 

3.  Determination  of  unit  flexible  airloads,  unit  incremental  flexible 
airloads,  and  CP's  on  each  strip,  and  wing  flexible-to-rigid  load 
ratios. 

4.  Determination  of  wing  flexible  lift  curve  slope. 

5.  Determination  of  wing  unit  flexible  shear,  bending  moment,  and 
torsion  along  the  reference  line  at  the  span  stations  specified  for 
weight  analysis. 


Wing  Unit  Rigid  Airload  Due  to  Angle  of  Attack 

Wing  unit  rigid  loadings,  rigid  lift  curve  slope,  and  CP's  are 
determined  as  follows. 

The  wing  platform  is  divided  into  10  equal -span  chordwise  strips 
between  the  side  of  the  body  station,  t^,  and  wingtip  station,  n * 1.0,  and 
one  strip  between  ngg  and  the  airplane  centerline  station,  n - 0,  as  shown 
in  Figure  3. 

For  subsonic  speeds  (M<1.0),  calculate  compressibility  factor,  B: 

B - (1  - mV/2  (1) 
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and  aspect  ratio  parameter,  BA/K: 
BA/K  = A • (B/K) 


(2) 


where  B/K  is  obtained  by  interpolation  of  B/K  versus  M data  in  Figure  4 
of  Volume  III,  Section  II. 

Calculate  compressible  sweep  parameter,  Afi: 


■ Arctan 


1 

tan  A.  „ - t- 

(}  • Ml] 

B 

LE  A 

\ 1 + X.  /J 

(3) 


Using  parameters  Ag,  BA/K,  and  interpolate  C f C/CL  C^y  data  in 
Table  1 of  Volume  III,  Section  II,  to  obtain  initial  C^C/Cg  C^y  values  at 
the  selected  rj  stations,  and  integrate  to  obtain  initial  unit  load  on  each 
strip;  i.e., 


n(IWBp/Cj  c 


r r I 

ll  lav, 


'(I)  INBD 


(4) 


Suimate  strip  loads  to  obtain  initial  total  wingload,  P^. 


11 

p = V ap  m 
wo  oKLJ 

i * i 


(5) 


Arbitrarily  change  carry-over  strip  No.  11  load,  APq  (11),  to  one-half 
initial  value  to  allow  for  lift  loss  due  to  body  interference;  i.e., 


APX  (11)  « AP0  (ll)/2 


(6) 


Sum  strip  loads  to  obtain  a new  total  wing  load: 


10 

pw  ' *pi  on  * I *po(I) 

I = 1 


(7) 
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The  unit  rigid  strip  loads  are  then  obtained  by  normalizing  a total  wing 
rigid  unit  load  of  1 per  side: 

For  strip  11, 

■ (^)(*P1(U))  (8) 

and  for  strips  1 through  10, 

' *PR.(I)  ■(V)(*V1))  (9) 

Using  parameters  Ag,  BA/K,  and  \ , interpolate  BC^K  data  in 
Table  2 of  Volune  III,  Section  II,  to  obtain  K value. 

The  initial  lift  curve  slope,  is  then, 

CLeW0  ’ 57.3  (BC^K)/(B/K),  per  raiian  (10) 

and  the  wing  rigid  lift  curve  slope  is  then, 

CL*WR  " (p^)CLoW0’  P*1*  radian  (11" 

For  supersonic  speeds  (M  >1.0),  calculate  the  compressibility 
factors : 


B - 

(M2  - 1)1/2 

02) 

Bn  » 

BCotALE 

(13) 

BA  « 

B-A 

w 

(14) 

Using  BA,  Bn,  and  \ , interpolate  (2r/Vi»b)  in  Table  3 of 
Volune  III,  Section  II,  at  selected  n stations,  and  integrate  to  obtain 
initial  load  on  each  strip;  i.e., 
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'(I)  OUTBD 


■ I'  Vi'  (2r/Vob)  dn 

7,1  (I)  INDB 


Sumnate  strip  loads  to  obtain  initial  total  wingload,  PWQ, 


WO 


l "o")  (•  ^ ) 


I - 1 


(15) 


(16) 


Hie  initial  value  of  wing  lift  curve  slope  is  then: 


f = a*  P 
LoWO  WO 


(17) 


Arbitrarily  change  carry-over  strip  11  load  AP0(11)  to  one-half 
initial  value  to  allow  for  lift  loss  due  to  body  interference;  i.e., 


APjCll)  - AP0(ll)/2 


(18) 


Sum  strip  loads  to  obtain  a new  total  wingload: 

10 

Pm  - APjfll)  • £ AP0(I) 

I = 1 


(19) 


Wing  unit  rigid  strip  loads  are  then  obtained  by  normalizing  to 
give  a total  wing  rigid  invt  load  of  1 per  side, 


*PR0fll)‘  (^f)(APl,n)) 


(20) 


and  for  strips  1 through  10 


apr,'i! 


APom| 


(21) 
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Wing  rigid  lift  curve  slope  is  then: 


S^wr  "(pwo)CLaWO  ^ 

Hie  column  matrix  of  unit  rigid  strip  loads  due  to  angle  of  attack, 

"•G-W  ' ApMi)  '23) 


where 


I • 1 through  10 
and  APRo  is  from  equation  9 or  21. 

CP's  of  unit  rigid  strip  loads  due  to  angle  of  attack  are  obtained 
as  follows: 

* For  strips  1 through  10,  the  CP  is  at  an  X/C  value  based  on  the 
mean  chord  of  the  strip.  The  (X/C)G  value  is  obtained  from  the 
(X/C) wing  versus  mach  number  plot  in  Figure  S of  Volume  III, 

Section  II. 

• For  strip  11,  the  CP  is  assumed  to  be  at  the  same  fuselage  station 
as  the  CP  of  strip  1. 

The  colunn  matrix  of  the  CP's,  AFWCP  (1,1),  for  strip  loads  due  to 
angle  of  attack  is  as  follows: 

AFWCP  (1,1)  - (X/C)  a (I)  (24) 


where 


1*1  through  10 


Wing  Unit  Rigid  Airload  Due  to  Wing  Flap  Deflection 

Wing  unit  rigid  strip  loadings  and  CP's  are  determined  as  follows. 
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Calculate  wing  span  ratios  for  outboard  and  inboard  ends  of  flap, 

bpo/b  * Yro/b/2  (25) 

bFI/b  " YFI/b/2  (26) 

Using  parameters  bj^/b  and  bpj/b,  interpolate  C£C/Cl0C^V  data  in 
Table  3 of  Volume  III,  Section  II,  at  selected  n values  for  bj^/b  and 
bpj/b  flap  span  ratios. 

Calculate  initial  span  loading  parameter  due  to  flap  deflection 
at  selected  n stations, 

^C,C/CUCAV)  ■ (C,C/CUCAV)  -|C|C/CUCAV)„ 

DW/b  °FI/b 

and  integrate  to  obtain  initial  load  of  each  of  the  strips. 


*PFO 


l 


n(I)  OUTBD 


<*C,C/CL„CAV>  *■ 


1(1)  INBD 

Simulate  strip  loads  to  obtain  flap  span  normalizing  parameter. 


(28) 


11 


^F  ^ APF0 
I = 1 


(29) 


Unit  rigid  strip  loads  due  to  flap  deflection  are  then, 


APRF(I)  “ Kgp  ’ APP0(I') 


(30) 


The  column  matrix  of  unit  rigid  strip  loads  for  the  lift  due  to 
flap  deflection,  AFW  (1,2),  is  as  follows: 


AFW(I,2)  = APrf(I) 


(31) 
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where 


I ■ 1 through  10 

CP's  of  unit  rigid  strip  loads  due  to  flap  deflection  are  obtained  as 
follows : 

* For  strips  1 through  10,  the  CP  is  at  an  X/C  value  based  on  the  mean 
chord  of  the  strip.  The  (X/C)p  value  is  obtained  from  the  (X/C)p 
versus  (Cf/Cy,)  data  in  Figure  6 of  Volume  III,  Section  II. 

* For  strip  11,  the  CP  is  assisted  to  be  at  the  same  fuselage  station 
as  for  strip  1. 

The  colirm  matrix  of  CP's  for  unit  rigid  strip  loads  due  to  flap 
deflection,  AFWCP  (1,2),  is  then: 

AF*CP  (1,2)  - (X/C)p  (I)  (32) 


where 


1*1  through  10 


Wing  Unit  Rigid  Load  Due  to  Vertical  Acceleration 


Uhit  rigid  strip  loads  due  to  vertical  acceleration  consist  of  the 
weight  of  the  strip  normalized  to  an  exposed  wing  load  of  -1  pound  per 
side;  i.e. , 


A?  (I) 
nzv  J 


Ws  (I) 


(33) 


where 


and 


Wg(I)  * strip  weight 
10 

Nor  ■ I ws  to 

I - 1 


(34) 


34. 


The  colunn  matrix  of  the  rigid  strip  loads  due  vertical  acceleration, 

AFW  (1,3)  is  then: 

AFW  (1,3)  - (I)  (35) 


where 


I ■ 1 through  10 


CP's  of  rigid  strip  loads,  (X/CW,  are  based  on  the  mean  chord  of  the 
strip,  as  shown  in  Figure  2.  The  column  matrix  of  the  CP's,  AFWCP  (1,3), 
is  then: 


AFVCP  (1,3)  - (X/C)w  (I) 


(36) 


where 


I ■ 1 through  10 


Wing  Flexible  Airload  Distributions 


The  methods  used  to  calculate  the  redistributed  wing  airloads  caused 
by  the  wing  torsional  and  bending  deflections  are  based  on  strip  theory. 
Redistributed  wing  loads  include  the  effects  of  the  wing  deflections 
resulting  from  (1)  lift  due  to  angle  of  attack,  (2)  lift  due  to  flap 
deflection,  and  (3)  inertia  load  due  to  vertical  acceleration. 
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Wing  Structural  Influence  Coefficients 


For  the  static  aeroelastic  analysis,  the  exposed  semispan  of  the  wing 
is  divided  into  10  evenly  spaced  strips,  as  shown  in  Figure  3.  Two  struc- 
tural influence  coefficient  points  are  placed  on  the  centerline  of  each 
strip,  one  at  X/C  forward  and  one  at  X/C  aft.  The  values  of  fuselage 
stations  (FS)  and  the  butt  line  (BL)  for  the  points  on  strip  I are  formed 
as  follows: 


BL  - BLgs  ♦ 


(37) 


" ^LE  ROOT  * BL  ' *“ALE 


(38) 


C “ CROOT  f1  ‘ b/2  (1 

»] 

(39) 

fsfwd  “ 

.W«FKD  • C. 

♦*IE 

(40) 

FSAFT  * 

Waft  ' c 

+ XLE 

(41) 

where 

CX/Q^d  - 0.15 
and 

WQm  - 0.65. 


The  coordinates  are  then  converted  into  the  swept  elastic  axis  Systran 
(Figure  3). 

x - (FS  - FS^q)  cos  A - BL  • sinA^  (42) 

y * (FS  - FS^q)  sin  ♦ BL  • cos  A^  (43) 

NOTE  The  elastic  axis  does  not  have  to  be  a constant-percent  chord  line. 
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Elements  of  the  SIC  matrix  are  formed  using  the  following  equation: 


SIC  (x,  y;  {,  n)  - f i\  ♦ 

fy 

xtJ0  cjd  * 


(44) 


Each  element  of  the  SIC  matrix  represents  the  deflection  at:  x,  y due 
to  a 1 -pound  vertical  load  at  f,  rf.  The  preceding  equation  is  integrated 
to  form  the  upper  right  triangle  of  the  SIC  matrix.  In  this  process  of 
integration,  the  value  of  y sometimes  exceeds  If.  In  that  case,  the  upper 
limit  of  integration  is  changed  to  W-  The  other  half  of  the  SIC  matrix  is 
formed  by  symmetry. 

SICi;j  - SICji  (45) 


For  the  static  aeroelastic  analysis,  a matrix  of  streamwise  slopes, 
SIC,  is  required.  This  matrix  is  formed  by  premulaplying  the  SIC  matrix 
by  a differential  operator  matrix  called  Dp  , which  is  formed  from  the 
structural  influence  point  geometry.  The  elements  of  D0  are  formed  as 
follows,  where  I goes  from  1 through  10. 


Dfl  (1,21-1)  - 


FSAFT^  ' FSFWD(I) 


(1,21)  - 


-1 


fsaftci>  * FWX> 

The  SIC  matrix  is  then  formed  as: 


(46) 

(47) 


[sic]  - [d„]  [sic] 
10,20  10,20  20,20 


(48) 
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Wing  Aerolastic  Distributions 

General  Method.  For  each  load  effect,  the  total  flexible  wing  load 
distribution  is: 


Where  the  incremental  aeroelastic  load  distribution  is  related  to  the 
flexible  load: 


{PAFLEx}“  q[A]  [^]  { PFLEX  } (5° 

In  this  simplified  solution,  the  aerodynamic  load  matrix,  A,  is  defined 
as  a diagonal  matrix  formed  from  the  a distribution  load  on  each  chordwise 
strip. 


[A]  - [a  LOAD] 


The  equation  for  the  aeroelastic  solution  is  then  formed  as: 


JPFLEX 


RIGID 


+ q 


[aLOADj 


Transpose  and  factor  out  PpLEX: 


f 1 

/ 

| I J - q aLQAD  SlC 

1 p 1 

| FLEX  j 

* ( PRIGID 

(51) 


(52) 


(53) 
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Let 


q 


a LOAD 


then, 


PRIGID 


(54) 


(55) 


Formulation  for  Flexible  Airloads  Program.  The  rigid  loads  distribution, 
AFW,  for  a , , an  WEIGHT,  are  distributed  to  the  SIC  points  by 

simply  beaming  the  load  on  each  strip,  using  the  corresponding  chordwise 
centers  of  pressure,  AFWCP. 


AFWCP(I,J)  - (X/C) 

R^ADS(n-l,J)  . AFW(I,J)  (x/c)fwd  , (x/c-)aft 


(56) 


RLQADS(2I,  J) 


-AFWCPCI  ,J)  ♦ (X/C)™ 
AFW(I’J)  (X/C)FWD  • WC)AFT 


(57) 


where 


and 


I ■ 1 through  10 

J * 1 a load 
J * 2 5F  load 
J s 3 weight 


Define  C as: 


C 


(58) 
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The  D-matrix  is  then  formed  a column  at  a time,  using  column  1 of  the 
RLQADS  matrix;  i.e.,  the  a load  distribution  and  the  3TC  matrix.  The  column 
K goes  1 through  20. 

D(2I-1,K)  - -C  • RLQADS (2 1 -1,1)  * SICRAR(I,K)  (59) 

D(2I,K)  - -C  * RLQADS(2I,1)  • SICRAR(I,K)  (60) 


where 


1*1  through  10 

The  D-matrix  is  completed  by  adding  1.0  to  each  diagonal  element. 
D(K,K)  - 1.0  + D(K,K) 

The  static  aeroelastic  equation  to  be  solved  is  then: 


T (FLEXIBLE)  ] f (RIGID) 1 
J [ FLOADS  J ” [RL0ADS  J 
20,20  20,3  20,3 


(61) 


This  equation  is  solved  for  the  flexible  loads  a column  at  a time,  using 
a general  least -squares  technique.  This  technique  is  used  because  it  is  very 
fast  and  will  work  with  even  a poorly  conditioned  matrix.  This  method  gives 
the  solution  to  a system  of  linear  equations  Bx  * C,  where  B is  an  N x M 
matrix  with  N > M,  and  C is  a column  vector  of  dimension  N.  In  the  over- 
determined case,  a least-squares  solution  is  obtained. 


The  method  finds  an  N x N matrix  R such  that  R is  orthogonal  (i.e., 
RTR*I)  and 


R(B.C) 


where  T is  an  M x M triangular  matrix  and  d is  an  M-dimensional  column  vector. 
The  number  a is  the  square  root  of  the  sum  of  the  squares  of  the  residuals. 
Solution  of  the  problem  is  then  the  solution  of  Tx*d,  and  is  found  by  back 
substitution. 
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Once  flexible  loads  at  the  SIC  points  have  been  determined,  flexible 
loads  on  the  strips,  AFWF,  are  formed  by  a simple  simulation;  the  delta 
flexible  loads,  AFWDF,  are  formed  by  subtracting  rigid  loads  from  flexible 
loads. 


AFWF(T,J)  - 

F LOADS  (21 -1,J)  + FLQADS(2I,J) 

(62) 

AFWDF  (I  ,J) 

- AFWF(I,J)  - AFW(I,J) 

(63) 

f 

where 

! 

| 

I - 1 through  10 

The  total  sunned  rigid  loads,  flexible  loads,  delta  flexible  loads, 
flexible-to-rigid  load  ratios  are: 

and 

10 

i 

Rigid 

F(J,1)  - ][  AFW(I,J) 

(64) 

t 

1-1 

10 

Flexible 

F(J,2)  - £ AFWF (I, J) 

(65) 

M 

A Flexible 

F(J ,3)  - F(J,2)  - F(J ,1) 

(66) 

F 

R 

p(j'4>  - w 

(67) 

Because  the  rigid  load  distribution  on  the  strips  is  used  for  delta 
flexible  loads,  the  chordwise  CP  on  each  load  strip  is  the  same  for  the 
delta  flexible  loads  as  it  is  for  the  rigid  a load. 
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The  wing  flexible  lift  curve  slope  is  determined  as  follows: 


CL  aWF 


CL  aWR 


10 

Apa  (11) 

♦I 

AFWF 

(i.i) 

1-1 

10 

Apa  (11) 

♦I 

AFW 

(1,1) 

1-1 

(68) 


where  APa  (11)  is  from  equation  8 or  20,  AFWF  (1,1)  is  the  column  matrix 
of  flexible  strip  loads  due  to  angle  of  attack,  and  AFW(I,1)  is  the  column 
matrix  of  rigid  strip  loads  due  to  angle  of  attack  from  equation  23. 


Wing  Unit  Shears  and  Moments 


Wing  unit  airload  shears  and  moments  for  the  flexible  wing  are  developed 
using  the  column  matrices  of  strip  loads  and  CP's.  The  average  unit  running 
normal  loading  on  a strip  is  the  strip  load  (AF  (I)  divided  by  the  strip 
width,  and  the  running  torsional  moment  is  the  running  normal  loading  times 
the  normal  distance  from  strip  CP  to  the  selected  load  reference  line  for 
weight  analysis.  The  average  running  normal  loading  and  running  torsional 
moment  are  applied  at  tj  station  for  the  midspan  at  the  strip. 

Running  loadings  on  the  exposed  wing  along  the  load  reference  line  are 
as  follows  (where  V (i)  is  the  rj  station  at  the  midspan  of  the  strip  I, 

UL(I)  is  the  unit  running  normal  loading  at  station  rj  (I) , and  UT(I)  is 
the  unit  running  torsional  moment  at  station  *?(!)): 


V(V  * figs  + 0.05  (1  - tj nc)  + 0.1  (I  - 1)  (1  -Tjnc) 


B S' 


BS 


-0.05  ♦ 1.05  BS  + (1  -tjbs)  1/10 


UL(I) 


AF  (I) 

0-1 


UT(I)  - -AXa(I)  • UL(I) 


(69) 

(70) 

(71) 
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The  term  AXA(I)  is  the  normal  distance  from  the  load  reference  line  to 
the  local  CP  and  is  determined,  using  the  methods  of  Volume  III,  Section  II, 
as  follows: 


AXaCI)  . CA(I)  [cXA/CA)CT  - CXA/CA)R]  (72) 

where 

[1  +»( I)  (1  -X)]  Cp  cosA. 

CA  (I) , - L (73) 

l-KA[l  - (X/C)R-  oca/ca)r] 

Ka  .1  (*i»ARCOSAR)  <74> 

(XA/CA)R  - (X/C)R  [l  - XA  (l  - (X/C)R)1  (75) 


(XA/CA)cp  ■ (X/C)^  [l  • XA  (l  - (X/C)cp)J  (76) 

Running  unit  normal  loading  inboard  of  the  side  of  the  body  (»j<  »»BS)  Is 
assumed  to  be  constant  and  is: 


UL„ 


AP  (ID 
’BS 


(77) 


where  AP(ll)  is  the  unit  airload  on  strip  11. 

Unit  running  torsional  moment  inboard  of  the  body  side  about  the  load 
reference  line  at  station  ij  is: 


- - AXa„  • UL„ 


where  AXAn  is  determined,  as  follows. 


(78) 
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The  CP  is  assumed  to  be  at  a constant  fuselage  station  equal  to  the 
fuselage  station  of  the  CP  at  strip  1.  Then,  when  n < ijgg, 

*xAn  " AXAn  (I  ^ + (n(i-i)'n)  (l)(sinAR)  (79) 

where  ^ is  from  equation  69 

and  AX.  is  from  equation  72 

" « ■ i) 

Unit  running  normal  and  torsional  loadings  for  the  specific  type  of 
distributions  are  obtained  using  the  following  strip  load  data: 

1.  For  flexible  lift  due  to  angle-of -attack  distribution,  use  the 
following : 

a.  In  equation  (70):  AF(I)  ■ AFWF(I,1)  from  equation  62 

b.  In  equation  (76):  (X/C)q,  ■ AFWCP(I,1)  from  equation  24 

c.  In  equation  (77):  AP(ll)  «APaR(ll)  from  equation  8 or  20, 

depending  on  mach  limber 

2.  For  rigid  lift  due  to  flap  deflection  distribution,  use  the 
following : 

a.  In  equation  (70):  AF(I)  * AFW(I,2)  from  equation  31 

b.  In  equation  (76):  (X/C)^,  ■ AFVCP(I,2)  from  equation  32 

c.  In  equation  (77):  AP(ll)  ■ APpR(ll)  from  equation  30 

3.  For  incremental  flexible  lift  due  to  flap  deflection  distribution, 
use  the  following: 

a.  In  equation  (70):  AF(I)  * AFWDF(I,2)  from  equation  63 

b.  In  equation  (76):  (X/C)^  » AFWCP(I,1)  from  equation  24 

c.  In  equation  (77):  AP(ll)  ■ 0 
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4.  For  incremental  flexible  lift  due  to  vertical  acceleration 
distribution,  use  the  following: 


a.  In  equation  (70):  AF(I)  ■ AFWDF(I,3)  from  equation  63 

b.  In  equation  (76):  (X/C)^p  ■ AFWCP(I,1)  from  equation  24 


c.  In  equation  (77):  AP(ll)  - 0 


Wing  unit  shear  and  moments  are  calculated  along  the  selected  load 
reference  line  at  the  n stations  selected  for  weight  analysis  and  are 
normalized  for  a surface  load  of  1 pound  per  side. 


Unit  span  normal  loading,  UL,  from  equations  70  and  77,  is  integrated 
to  obtain  the  unit  shear  values,  USZN,  at  the  selected  ij  stations;  i.e., 


USZN  - Kjj2 


1.0 

/ (ULN)d„ 


(80) 


where 


*uz- 


1.0 

TO 

J (ULN)dij 

0 


(81) 


Unit  shear,  USZ,  is  integrated  to  obtain  the  unit  bending  moment, 


UMXTN: 


UMXTN 


b 

2 cos  Ar 


(USZN)d 


(82) 
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Unit  running  torsional  moment,  UT,  from  equations  71  and  78,  is 
integrated  to  obtain  the  unit  torsional  moment,  UMYTN: 

1.0 

UM YTN  - J mi)drf  (83) 

n 

In  equations  80  through  83,  the  N in  ULN,  USZN,  UMXTN,  UMYTN,  and  UTN 
is  used  to  designate  the  type  of  distribution;  i.e., 

N - 1 for  flexible  lift  due  to  angle  of  attack 

N * 2 for  rigid  lift  due  to  flap  deflection 

N ■ 3 for  incremental  flexible  lift  due  to  flap  deflection 

N « 4 for  flexible  lift  due  to  vertical  acceleration 

Wing  unit  airloads  and  CP's  for  each  of  the  four  types  of  airload 
distribution  are  determined  as  follows. 

Total  wing  unit  airload  per  side,  UPZWN;  wing  unit  pitching  moment  per 
side,  UMYWN,  at  the  intersection  of  the  wingload  reference  line  and  the  A/P 
centerline;  and  CP  of  the  wing  unit  airload,  DXWN,  measured  aft  from  the 


theoretical  wing  apex,  are  as  follows: 

UPZWN  - (USZN)  _ 0 - 1.0  (84) 

UMYWN  - (UMYTN) ^ m q cos  Ar  - (UMXTN) ,,  , Q sinAR  (85) 

DXWN  « CR(X/C)R  - (UMYWN) /(UPZWN)  (86) 

Exposed  wing  unit  loads  of  the  intersection  of  the  wingload  reference 
line  and  the  body  side  in  the  body  reference  system  are  as  follows: 

USZW(B)N  - (USZN)  „ gg  (87) 

UMXW(B)N  - (UvDCIN)  nfiS  cosAR  ♦ (UMYTN)  ^ sinAR  (88) 

UMYW(B)N  - (UMYTN)  n RS  cos  AR  - (UMXTN)  „ K sinAR  (89) 
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NOTE  Also 


USZW(B)N 


(USZN)  g 

cuszn)  ; “ 


UPZWN 


(90) 


mm  - C^0t/C)R  * Y^tanA,  - (91) 

W(«N  ' YBS  * SI  (92) 

Wing  unit  carry-over  load  on  the  body,  UPZB(W)N,  and  its  CP,  DXB(W)N, 
are  determined  as  follows: 

UPZB(W)N  » (USZN)  „ mQ  - (USZN)  „ ^ (93) 


UMYB(W)N  - ftJMVTN  „ _Q  - LMYlNiijs]  C0SAR  - 


iwmi 


, -o  - 1*®CTO ■> BS  - (USZN->BS) 


\ cos  ar  / 


sinAR  (94) 


DXB(W)N 


CR(X/C) 


UMYB(W)N 
R ‘ UPZltwjN 


(95) 


EMPENNAGE  UNIT  AIRLOAD  DISTRIBUTIONS 

The  unit  airload  distributions  on  the  horizontal  and  vertical  tail 
surfaces  are  determined  using  the  methods  described  in  Volume  III, 
Section  II. 
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MT  LIMIT  AIRLOADS 


Gross  limit  airloads  on  airplane  components  and  airplane  inertia  factors 
are  determined  for  specific  types  of  flight  conditions,  using  methods  similar 
to  those  described  in  Voltne  III,  Section  II.  Revisions  are  made  to  include 
the  effects  of  wing  flexibility,  consisting  of  the  following: 

1.  Wing  airload  is  divided  into  four  parts:  (1)  flexible  load  due  to 

angle  of  attack,  (2)  rigid  load  due  to  flap  deflection,  (3)  incre- 
mental flexible  load  due  to  flap  deflection,  and  (4)  incremental 
flexible  load  due  to  vertical  acceleration. 

2.  Airplane  balance  equations  are  expanded  to  incorporate  added  wingload 
parts. 

3.  Flexible  wing  lift  curve  slope  is  used  to  detexmine  wing  incremental 
gust  load. 


BALANCED  MANEUVER  CONDITION 

The  balanced  maneuver  condition  is  a flight  condition  where  the  aircraft 
is  trimmed  (balanced)  by  a horizontal  tail  or  canard  load.  Incremental  loads 
due  to  flap  deflection  are  determined  only  for  low  subsonic  conditions. 

Wing  airloads  consist  of  the  following: 

PZW1  ■ Flexible  lift  due  to  angle  of  attack 

PZW2  ■ Rigid  lift  due  to  flap  deflection 

PZW3  * Incremental  flexible  lift  due  to  flap  deflection 

PZW4  ■ Flexible  lift  due  to  vertical  acceleration 

Their  corresponding  CP's  in  terms  of  fuselage  station,  are: 


XW1  « + DXW1 

(96) 

XW2  - + DXW2 

(97) 

XW3  - xRW  + DXW3 

(98) 

XW4  ■ + DXW4 

(99) 
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DXWN  values  are  obtained  from  equation  (86)  for  each  type  of  distribution. 

For  a given  mach  No.,  dynamic  pressure,  load  factor  and  flap  deflection, 
wingloads  are  determined  as  follows: 


PZW1  is  obtained  from  the  balance  equation  (110) 
PZW2  " 5773  KCF  CL»WRqSW 


(100) 


Where 


6^  ■ flap  deflection  in  degrees 


- 


CF  dC./do 

and  is  obtained  by  interpolation  of  the 

Kq.  versus  (Cp/C^)  data  in  Figure  6 of  Volume  III,  Section  II 

Kgp  is  obtained  from  equation  29 

CLaWR  is  obtained  *Tom  equation  11 

q ■ dynamic  pressure,  psf 
Sw  * theoretical  wing  area 

PZW3  - (PZW2)  (USZW2)s0b  - l.oj 


(101) 


Where  (USZW2)S0B  is  the  unit  rigid  wingload  at  the  body  side  due  to  flap 
deflection  and  is  obtained  from  equation  87. 

F/R  ■ matrix  element  F(2,4)  from  equation  67 


PZW4  - 2.0  (Nz)  (Ww0p)  1.0 


(102) 
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Where  is  the  specified  load  factor  and 


WW0P  " exP°se<*  wing  weight  per  side 


-jr  - matrix  element  F(3,4)  from  equation  67 


Body  nose  load,  PrsJ,  is  as  defined  as  follows: 


PZN  * ft2wRN  q/l44 


0.043633  a q 


(103) 


and 


Where  % is  the  maximum  nose  radius  (inches) , a is  the  angle  of  attack 


PZWl/qSyy 

CLoWF 


(104) 


Where  Cl^^F  *s  t^ie  flexihle  wing  lift  curve  slope  and  is  obtained  from 
equation  68. 

Substituting  equation  104  in  103: 


ZN 


0. 043633 (PZW1)  R^ 
SW  CLoWF 


(105) 


CP  of  the  body  nose  load  is: 


XN  - X 


* (‘<*--*-2 

\ *v 


(106) 


Horizontal  tail  load  is  P m and  is  obtained  from  the  balance  equa- 
tion 111.  CP  of  the  horizontal  tail  is: 


XH  = X + DXH 
HH 


(107) 
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Where  Xrh  is  the  fuselage  station  of  the  leading  edge  of  the  horizontal 
tail  theoretical  root  chord,  DXH  is  the  X-distance  aft  of  the  root  chord  lead* 
ing  edge  to  the  CP  and  is  obtained  from  equation  50  in  Volume  III,  Section  II. 

Airplane  balance  equations,  where  Q ■ 0,  are  us  follows: 


EMy  - 0 


(PZW1) 


(XCG  - XW1)  + (XCG  - XN) 


(0.043633  Rjj2  \ 
^ CLoWF  / 


PZW2  (XCG  - XW2)  + PZW3  (XCG  - XW3)  + 
PZW4  (XCG  - XW4)  + PZH  (XCG  - XH)  * 0 


(108) 


(109) 


Solving  equations  108  and  109  for  PZW1  and  PZH, 

(NZW)F  - (PZW2)  (C*F)  - (PZW3)  (D-F)  - (PZW4)  (E-F) 


PZW1  » 


B - AF 


PZH  - NZW  - A(PZW1)  - PZW2  - PZW3  - PZW4 


(110) 

(111) 


where 


A - 1 + 


0.043633  R^2 
CLoWF 


B - (XCG  - XW1)  + (XCG  - XN) 


(112) 


(113) 


c - 

(XCG 

- XW2) 

(114) 

D - 

(XCG 

- XW3) 

(115) 

E - 

(XCG 

- XW4) 

(116) 

F - 

(XCG 

- XH) 

(117) 

PZN 

- (A- 

1)  PZW1 

(118) 

PITCHING  ACCELERATION  CONDITION 


Hie  pitching  acceleration  condition  is  an  arbitrary  condition  where  a 
specified  value  of  pitching  acceleration  is  caused  by  an  incremental  hori- 
zontal tailload  and  is  superimposed  on  a balanced  maneuver  condition  such  that 
the  resulting  load  factor  is  one-half  the  design  limit  positive  maneuver  load 
factor. 


The  incremental  horizontal  tailload  required  to  produce  the  specified 
pitching  acceleration  0.  is: 


APZH  - -12 


Q Iy 

XH  - XCG 


(119) 


Where,  Iy  ■ airplane  pitching  moment  of  inertia,  slug  ft  , and  XH  is 
obtained  from  equation  107. 


The  maneuver  load  factor,  N_„  is  determined  as  follows: 

ZM 


..  nzl  apzh 

M S - ■ ...  ■■ 

ZM  2 W 


(120) 


Where  N_.  is  the  design  limit  maneuver  load  factor 


W * airplane  weight 

Equations  96  through  118,  with  the  following  exceptions,  are  then  used  to 
determine  the  component  loads: 


(121) 
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N-,  F - (PZW4)  (E-F) 

PZW1-  B - AT j 

1 

1 

(122) 

PZH  - N^  - A(PZW1)  - PZW4  + APZH 

(123) 

VERTICAL  GUST  CONDITION 


The  vertical  gust  condition  consists  of  a ±50  fps  vertical  gust  encounter 
superimposed  on  a l.Og  trained  condition. 

Equivalent  airspeed  (Vg) , airplane  mass  ratio  (hO  , and  the  gust  allevia- 
tion factor  (Kg)  are  determined  using  equations  115  through  119 of  Volume  III, 
Section  II,  at  the  specified  mach  number  and  altitude. 


Component  loads  PZWlNz  ■ i.o  and  PZ%-  « , «,  for  the  l.Og  trimmed  con- 
dition are  determined  using  equations  (110J  ana  (111),  where  Nz  * 1.0  and 

6f  - 0. 

Airplane  normal  load  factor,  Nz,  is  determined  as  follows: 


NZW 


l™1 N 


+ APZW1  \ - 2.0  N„ 


Z - 1.0 


^WOP^ 


(H 


+ APZH 

1.0 


(124) 


or: 


N„* 


A(PZW1N  + APZW1)  + PZ^ 


+ APZH 


1.0 


w*2'0  (I'1) 


(125) 


where  PZWIn^  « i q ^ obtained  using  equation  (110) , PZHnz  b obtained 
using  equation  (ill) , A is  obtained  using  equation  (112) , F/R  is  the  same  as 
that  used  in  equation  (102),  and  incremental  gust  loads  are: 


4PZW1  - 0.100354  Kg  Cuwp  ^ V£ 


(126) 


APZH-  0.100354  Kg 


S3 


Where  Cl^  is  obtained  using  equation  (68)  and  Cl^  is  obtained  using 
equation  19  or  23  of  Volume  III,  Section  II. 

Limit  airload  on  the  components  are  then  determined  as  follows: 


PZW1  - PZWL.  + APZW1 

^Z  - 1.0 

(127) 

PZW4  - - 2.0  Nz  (W^I  - l) 

(128) 

PZH  - PZH..  + A PZH 

^Z  - 1.0 

(129) 

PZN  - (A  - 1)  PZW1 

(130) 

Airplane  pitching  acceleration,  Q,  is  determined  as  follows: 

Q = - [(XN-XCG)  PZN  + (XW1  - XCG)  PZW1  + (XW4  - XCG)  PZW4  + 

(XH  - XCG)  PZH]  / 12  Iy  (131) 

LATERAL  GUST  CONDITION 

The  lateral  gust  condition  consist  of  a 50  fps  lateral  gust  encounter 
superimposed  on  a 1.0  g trimmed  condition. 

Component  limit  airloads  PZW1,  PZW4,  PZH,  and  PZN  are  determined  using 
equations  110,  102,  111,  and  118,  respectively,  with  Nz  = 1.0  and  6p  = 0. 

Component  limit  airloads  Pyjv  and  PyvT  are  determined  using  equations  138 
and  139,  respectively,  from  Volume  III,  Section  II. 

Airplane  yawing  acceleration,  R,  and  side  load  factor,  Ny,  are  deter- 
mined using  equations  140  and  141 , respectively,  from  Volume  III,  Section  II. 


YAWING  ACCELERATION  CONDITION 

The  yawing  acceleration  condition  is  an  arbitrary  condition  where  a 
specified  value  of  yawing  acceleration  is  caused  by  a load  on  the  vertical 
tail  and  is  superimposed  on  a 1.0  g trinmed  condition. 


54 


Component  limit  airloads  PZWl,  PZW4,  PZH,  and  PZN  are  determined  using 
equations  110,  102,  111,  and  118,  respectively,  with  ■ 1.0  and  6p  ■ 0. 

Equations  142  and  143,  of  Volume  III,  Section  II,  are  used  to  determine 
vertical  tail  limit  airload,  P and  airplane  side  load  factor,  Ny, 
respectively.  ^ 


LIFTING  SURFACE  LIMIT  AIRLOAD  SHEARS  AND  MOMENTS 


Lifting  surface  unit  airload  shears  and  moments  have  been  determined  (for 
each  type  of  distribution)  such  that  the  panel  load  (per  side  including  the 
carry-over  load  on  the  body)  is  equal  to  unity.  Therefore,  limit  airload 
shears  and  moments  are  obtained  by  multiplying  the  unit  values  by  the  limit 
panel  loads. 


WING  LIMIT  AIRLOAD  SHEARS  AND  MCMENTS 

Wing  airload  shears,  bending  moments,  and  torsional  moments  at  the 
selected  n stations  for  weight  analysis  along  the  load  reference  line  are 
determined  for  a specific  flight  condition,  as  follows: 

Ya  - (b/2)/Cos  Ar  - Y/CosAr  (132) 


S * [PZWl(USZl)  + PZW2(USZ2)  + 

TV 

PZW3(USZ3)  + PZW4  (USZ4)]  +2  (133) 


My  *=  [PZWl  (UNDO!)  + PZW2  (UMXT2)  ♦ PZW3(UMH3)  + 

PZW4  (UMXT4)]  +2  (134) 


NL  - [PZWl  (UMYT1)  + PZWZ  (UMYT2)  + 

PZW3  (UMYT3)  + PZW4  (UMYT4)]  + 2 (135) 


where  the  unit  shears  and  the  unit  bending  and  torsional  moments  are  obtained 
from  equations  80,  82,  and  83,  respectively,  for  each  type  of  distribution. 
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Limit  airload  shear  and  moments  of  the  side  of  the  body  stations,  Ygg,  in 
the  body  reference  system  are  obtained  using  unit  values  from  equations  87,  88, 
and  89  for  each  type  of  distribution. 


[PZW1  (USZW(B)l)  + PZW2  (USZW(B)2)  ♦ 

(USZW(B)3)  + PZW4  (USZW(B)4)]  + 2 (136) 

[PZW1  (UMXW(B)l)  + PZW2  (IMXW(B)2)  + 

(UMXW(B)3)  + PZW4  (UMXW(B)4)  ] + 2 (137) 

[PZW1  (LMYW(B)l)  + PZW2  (UMYW(B)2)  + 

(LMYW(B)3)  + PZW4  (UMYW(B)4)]  + 2 (138) 

The  total  limit  carry  over  airload  on  the  body,  PZB(W) , and  its  CP, 

XBB(W) , are  determined  as  follows  using  unit  values  from  equations  93 
and  95: 

PZB(W)  « PZW1  (UPZB(W)1)  ♦ PZW2  (UPZB(W)2)  + 

PZW3  (UPZB  (W)  3)  + PZW4  (UPZB(W)4)  (139) 

XBBOV)  ■ + [PZW1  (UPZB(W)1)  DXB(W)1  ♦ 

PZW2  (UPZB (W)  2)  DXB  (W)  2 + PZW3  (UPZB (W)  3)  DXB(W)3  + 

PZW4  (UPZB(W)4)  DXB(W)4]  h-PZB(W)  (140) 

EMPENNAGE  LIMIT  AIRLOAD  SHEARS  AND  MOMENTS 

Limit  airload  shears  and  bending  and  torsional  moments  are  determined 
using  equations  153  through  182,  in  Volume  III,  Section  II. 


'(SOB) 

PZW3 


(SOB) 

PZW3 


(SOB) 

PZW3 
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Section  III 


COMPUTER  PROGRAM  DESCRIPTION 


GENERAL  DESCRIPTION 


The  flexible  airloads  stand-alone  program  BFCNTL,  has  been  developed  to 
determine  airplane  component  limit  airloads,  including  effects  of  wing 
flexibility  for  specific  flight  case  conditions.  The  limit  airloads  consist 
of  the  airload  and  center  of  pressure  (CP)  for  each  airplane  component  and 
airload  shear,  bending  moment,  and  torsion  distributions  on  wing  and 
empennage  surfaces. 

Punched  card  output  data  are  in  a format  that  is  compatible  for  use  as 
an  optional  external  input  to  the  SWEEP  program.  Input  data  are  generated  by 
the  SWEEP  data  generation  program  as  punched  card  data  such  that  the  operation 
of  the  stand-alone  program  requires  no  further  effort  other  than  setting  up 
the  program  decks  and  the  preparation  of  a control  data  card. 


PROGRAM  FUNCTIONS 


The  flexible  airloads  main  program,  BFCNTL,  utilizes  three  main  subrou- 
tines and  nine  other  subordinate  subroutines.  The  main  subroutines  are 
USPANF,  BNLDSF,  and  SPABMF. 

Subroutine  USPANF  is  used  to  determine  lifting  surface  unit  airload 
shears,  moments,  CP's  and  lift  curve  slopes  for  a specified  mach  nunber  and 
altitude.  For  the  empennage  (horizontal  and  vertical  tail)  surfaces,  the 
methods  used  are  described  in  Volume  III,  Section  II.  For  the  wing,  methods 
used  are  described  Section  II  herein. 

Subroutine  BNLDSF  is  used  to  determine  the  gross  limit  airload  and  CP  on 
each  of  the  airplane  components  and  the  airplane  inertia  factors  for  the 
specified  case  conditions.  Methods  employed  are  described  in  Section  II. 

Subroutine  SPABMF  is  used  to  deteimine  limit  airload  shear,  bending 
moment,  and  torsion  distributions  on  lifting  surfaces  for  a specific  flight 
condition.  For  empennage  surfaces,  the  methods  described  in  Volume  III, 
Section  II,  are  used.  For  the  wing,  methods  used  are  described  in 
Section  II  herein. 
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The  purposes  of  other  subordinate  subroutines  are  as  follows: 


• DECRD 

• RERDAT 

• C0DIM2 

• FC0EM2 

• ATM0S 

• WFLEX 

• FLXSIC 

• GLSQ 

• MATRIT 


Reads  and  stores  input  data  in  assigned  dimension  region 

Rearranges  and  stores  input  data  in  assigned  locations 

Curve  fitting  interpolation  routine  for  determination  of 
a value  on  a single  curve 

Curve- fitting  interpolation  routine  for  the  determination  of 
a value  from  a family  of  curves. 

Detexmines  the  standard  atmosphere  density  and  the  speed  of 
sound  for  a given  altitude 

Determines  the  flexible  wing  airload  distribution  data  for 
given  types  of  rigid  wing  and  load  data 

Deteimines  wing  structural  influence  coefficient  data 
required  by  WFLEX 

Gives  the  solution  to  a system  of  linear  equations. 

Routine  for  printing  matrices  by  rows  or  colunns 


The  calling-called  matrix  for  the  program,  showing  the  interdependence  of 
the  subroutines,  is  shown  in  Figure  4. 


MAIN  PROGRAM  - BFCNTL 


The  main  program,  BFCNTL,  performs  the  following  functions: 

1.  Reads  in  input  ND,  BC,  and  BF  data  arrays  of  which  BC  and  BF  are  the 
punched  card  outputs  from  the  SWEEP  II  data  generation  program. 

2.  Rearranges  BC  input  array  and  sets  up  flight  condition  data  for 
subroutine  usage. 

3.  Calls  subroutines  and  transfers  data  required  to  perform  the  airload 
calculations . 

4.  Punches  card  output  (a&id  prints)  data  containing  component  limit 
airloads  and  CP's  and  the  wing  and  empennage  limit  airload  shear 
and  moment  distributions. 
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Figure  4.  Calling-called  matrix  for  flexible  airloads  stand-alone  program. 
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Methods  described  in  Section  II  are  used  in  the  coding  of  the  program. 
The  logical  flow  chart  for  the  main  program,  BFCNTL,  is  shown  in  Figure  S, 
and  detail  flow  charts  and  program  listing  are  presented  in  Appendix  A. 


PROGRAM  DECK  SETUP 

The  program  deck  setup  is  illustrated  in  Figures  6 and  7.  The  total 
deck  setup  must  follow  the  blocked  order  shown.  Twelve  subordinate  subrou- 
tines may  be  arranged  in  any  order,  but  the  total  subroutine  block  must  be 
immediately  behind  the  main  program  deck.  When  only  one  case  is  to  be  run, 
the  execute  card  must  follow  the  last  card  of  DP  data  set. 

When  multiple  cases  are  to  be  run,  subsequent- case  data  must  be  arranged 
as  shown  in  Figure  7 and  placed  immediately  behind  the  first  case.  The 
execute  card  must  then  follow  the  last -case  DP  data  termination  card. 


INPUT  DATA 

Input  data  for  the  program  consist  of  a user-prepared  ND  array  and  the  BC 
and  BF  arrays  produced  by  SWEEP  data  generation  option. 

The  input  ND  data  array  is  an  integer  array  on  one  card  containing  con- 
trol factors  shown  in  Table  1.  The  flexible  loads  program  uses  some  coding 
procedures  which  are  identical  to  those  contained  in  the  airloads  estimation 
program  in  Volume  III.  Consequently,  some  of  the  control  factors  listed  in 
Table  1 are  identical  to  those  in  Volume  III.  To  simplify  the  coding  proce- 
dures, identical  control  factors  are  located  in  the  same  card  columns;  some 
columns  are  left  blank;  and  some  columns  contain  relocated  and  new  control 
data. 


Figure  1 and  Table  2 show  basic  conditions  produced  for  each  of  the 
various  case  numbers.  Any  basic  condition  not  specifically  defined  by  the 
chosen  case  nunber  will  not  be  produced,  even  if  the  type  factor  ND(28) 
through  ND(36))  is  entered  as  yes.  These  type  factors  are  used  for  selection 
or  rejection  of  only  those  basic  flight  conditions  described  by  the  chosen 
case  number. 

BC  data  array  contains  (1)  airplane  design  weights , centers  of  gravity 
locations,  and  moments  of  inertia,  (2)  design  limit  load  factors,  (3)  design 
speed- altitude  points,  (4)  airplane  dimensional  and  geometric  data,  and 
(5)  span  stations  on  the  lifting  surfaces  selected  for  weight  analysis. 

These  data  are  on  punched  cards  in  the  E- format  that  is  compatible  with 
decimal  read  subroutine  DECRD.  Input  BC  array  is  shown  in  Table  3. 
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Main  program  BFCNTL 


Figure  3.  Logical  flow  chart  for  flexible  airloads  stand-alone  program. 
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NO  factors  card 


Figure  7.  A subsequent -case  data  block  setup. 
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BF  data  array  contains  (1)  design  speeds  at  20,000  feet,  (2)  wing 
elastic  axis  location,  (3)  exposed  wing  strip  weights  and  center  of  gravity 
(X/C)  locations,  and  (4)  wing  spanwise  variation  of  the  stiffness  parameters 
El  and  GJ.  These  data  are  shown  in  Table  4,  locations  1 through  106,  and  are 
on  punched  cards  in  the  same  format  as  the  BC  data. 


SEMIPERMANENT  DATA 

The  program  deck  setup  contains  data  sets  DT,  DB,  DF,  and  DP,  which 
comprise  semipermanent  (fixed)  data.  These  data  sets  contain  aerodynamic  data 
used  by  subroutine  USPANF  in  the  determination  of  lifting  surface  unit  airload 
distributions.  These  data  sets  are  the  same  as  that  contained  in  the  SWEEP 
airloads  module  data  bank. 

DT  data  set  array  is  shown  in  Table  5.  This  data  set  contains  n-span 
stations  and  taper  ratios  applicable  to  spanwise  loading  parameters  of  DB  and 
Iff  data  sets,  and  contains  the  variation  with  mach  number  of  compressible 
section  lift-curve  slope  parameter,  B/K. , and  section  CP,  X/C,  for  the  wing 
and  empennage. 

DB  data  set  array  is  shown  in  Table  6.  This  data  set  contains  subsonic 
span  loading  parameter  (Cj  C/CL  C^y)  * variation  of  compressible  lifting  sur- 
face lift- curve  slope  parameter  (BCj^/K)  versus  compressible  sweep  parameter 
(Xg)  and  aspect  ratio  parameter  (BA/Ki)  for  16  combinations  of  span 
station  (n)  and  taper  ratio  (k). 

The  DF  data  set  array  is  shown  in  Table  7.  This  data  set  contains 
section  CP  (X/C)  and  flap  lift  effectiveness  parameter  (Kcf)  variation  with 
flap  chord  ratio  (Cp/Cw)  and  flap  lift  spanwise  loading  parameter 
(C  g C/Cl<>CAV)  variation  with  span  station  ( r\ ) for  10  flap  span  ratios 
(bp/tty) . 

DP  data  set  array  is  shown  in  Table  8.  This  data  set  contains 
supersonic  span  loading  parameter  (2r/Vcrb)  variation  with  sweep  parameter 
(Bm)  and  aspect  ratio  parameter  (BA)  for  16  combinations  of  span  station  ( n ) 
and  taper  ratio  ( k ) . 

Numeric  values  of  DT,  DB,  DF,  and  DP  data  contained  in  the  program  data 
bank  are  presented  in  Volume  III,  Section  II. 
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SUBROUTINE  DESCRIPTIONS 


SUBROUTINE  RERDAT 

Subroutine  RERDAT  reads  in  and  rearranges  the  input  BC  data  set 
described  in  Table  3 to  data  locations  shown  in  Table  9.  The  subroutine  also 
interpolates  El  and  GJ  data  BF(47)  through  BF(106)  to  obtain  El  and  GJ  data 
for  subroutine  WFLEX,  as  shown  in  data  locations  BF(114)  through  BF(173)  of 
Table  4.  BF  data  items  BF(107)  through  BF(113),  BF(174),  and  BF(175)  are 
added  to  the  input  data  list.  The  rearranged  BC  data  locations  and  the 
complete  BF  data  sets  are  then  used  in  the  main  program  and  the  subprogram. 

Subroutine  RERDYT  logical  flow  chart  and  its  program  listing  are  shown  in 
Appendix  A. 


SUBROUTINE  DECRD 
Entry  Name 


DECRD  provides  the  facility  for  reading  a variable  number  of  pieces  of 
real  data  from  the  input  device  and  storing  them  in  specified  elements 
(either  sequential  or  nonconsecutive)  or  an  array. 


Usage 


CALL  DECRD  (APR) 

APR  - The  name  of  the  real  array  to  be  read. 

This  routine  is  particularly  helpful  in  programs  in  which  the  number  of 
input  elements  varies  from  case  to  case.  Only  information  specified  is 
actually  read  into  storage;  remaining  elements  of  the  array  are  unchanged. 

Data  are  usually  written  on  the  form  Fortran  Fixed  10  Digit  Decimal 
Data.  Each  card  must  contain  an  index:  an  integer  written  in  columns  2 

through  12.  The  five  data  fields  of  12  columns  each  (columns  13  through 
72)  contain  input  data  of  the  real  type.  However,  any  data  field  may  be  left 
blank  to  indicate  that  the  corresponding  location  in  core  is  not  to  be 
changed.  Columns  73  through  BO  contain  the  identification. 

The  index  defines  the  location  of  the  first  piece  of  data  on  the  card 
within  the  array  specified  as  the  argument.  This  integer  must  be  written  to 
the  extreme  right  of  the  field.  If  the  name  of  the  array  is  not  subscripted 
in  the  CALL  statement,  the  index  can  be  considered  equivalent  to  the  subscript 
of  a one -dimensioned  array.  For  example,  if  the  argument  in  the  CALL  is  the 
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nonsubscripted  array  name,  ARh,  and  the  index  is  10,  the  first  piece  of  data 
on  the  card  (columns  13  through  24)  will  be  read  into  ARR  (10);  the  third 
piece  of  data  (columns  37  through  48)  will  be  read  into  ARR  (12). 

For  an  array  with  multiple  subscripts , the  index  should  be  computed  so 
that  the  particular  element  can  be  defined  by  a single  number.  The  index 
may  not  be  zero  or  blank. 

All  data  items  must  be  of  the  real  type;  they  are  written  following  the 
rules  for  input  with  E-type  format  specification.  If  an  exponent  is  written, 
it  must  be  at  the  extreme  right  of  the  field. 

1.  If  the  number  is  written  without  either  an  exponent  or  a decimal 
point,  the  point  is  assumed  to  be  at  the  extreme  right  of  the  field 
(as  if  read  with  an  E12.0  format). 

2.  If  the  decimal  point  is  explicitly  written,  the  number  may  be 
positioned  anywhere  in  the  field. 

3.  If  no  decimal  point  is  written  but  an  exponent  is  furnished,  the 
point  is  assumed  to  be  immediately  to  the  left  of  the  exponent. 

When  a field  is  left  blank,  no  information  is  read  into  the  location 
corresponding  to  this  field;  the  information  already  in  this  location  is 
unaltered.  A negative  zero  is  read  as  zero. 

Reading  is  terminated  by  putting  a negative  sign  in  column  1 of  the 
last  card  to  be  read. 


Error 


If  card  columns  2 to  12  are  blank  or  zero,  the  comment,  "DECRD  ER. 

CARD  ■ (bad  card  image)"  is  printed  and  execution  of  the  job  is  terminated. 
If  a field  contained  in  columns  13  through  72  cannot  be  converted  by  the 
specified  format,  an  error  code  is  printed  and  execution  of  the  job  is 
terminated. 


Method 

Data  fields  of  each  card  are  converted  twice,  using  two  formats, 
SF12.0  and  10A6.  The  A- type  conversion  is  used  to  check  for  blank  fields. 
If  the  field  is  not  blank,  the  result  of  the  E-type  conversion  is  stored 
in  the  proper  element  of  the  specified  array.  After  reading  each  card,  a 
test  is  made  for  a negative  sign  in  the  first  field;  reading  is  terminated 
if  the  sign  is  negative.  The  numerical  data  cards  processed  by  \i  s \b- 
routine  is  discussed  in  Section  I of  Volume  IX. 


66 


SUBROUTINE  USPANF 


Subroutine  USPANF  is  used  to  determine  the  lifting  surface  unit  airload 
shears,  moments,  CP's  and  lift -curve  slopes  for  a combination  of  mach 
number  and  altitude  specified  by  the  main  program. 

Methods  used  for  the  horizontal  and  vertical  tail  surface  are  the  same 
as  those  described  in  Volume  III,  Section  II,  for  the  SWEEP  program. 

Methods  used  to  develop  the  wing  unit  airload  data  include  the  airload 
distributions  caused  by  wing  bending  and  torsional  deflections  resulting 
from  (1)  lift  due  to  angle  of  attack,  (2)  lift  due  to  flap  deflection,  and 
(3)  inertia  load  due  to  vertical  acceleration.  The  applicable  methods  are 
described  in  Section  II  herein. 


FUNCTION  C0DIM2 

Function  C0DIM2  is  an  interpolation  routine  for  the  determination  of  a 
point  on  a single  curve  fitted  through  four  points. 


FUNCTION  FC0EM2 

Function  FO0EM2  is  an  interpolation  routine  for  the  determination  of  a 
point  from  a family  of  fitted  curves.  The  subroutine  utilized  C0DIM2  for  each 
family  curve  interpolation. 


SUBROUTINE  ATM0S 

Subroutine  ATM0S  is  used  to  determine  the  1962  U.S.  Standard  Atmosphere 
density  and  speed  of  sound  for  a given  geometric  altitude.  The  methods 
employed  can  be  followed  in  the  logical  flow  chart  and  the  program  listing 
shown  in  Appendix  A. 


SUBROUTINE  WFLEX 

Subroutine  WFLEX  solves  the  wing  static  aeroelastic  problem  for  the  a, 
6 FLAP,  and  vertical  inertia  load  effects.  The  subroutine  uses  distributed 
rigid  airload  data,  along  with  flight  condition  data  and  structural  influe- 
ence  coefficients  from  subroutine  FLXSIC.  The  logical  flow  chart  and  the 
program  listing  are  shown  in  Appendix  A. 
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Usage 


All  input -output  data,  except  printed  results  and  diagnostic  printed 
results,  are  passed  through  the  call  statement. 

CALL  WFLEX  (AFW,  AFWCP,  AFWF,  AFWDF,  F,  YEIGJ,  El,  GJ,  CR,  B02,  BLBS, 
ANGLE,  FSLERT,  ANGERA,  FSEA0,  WLAMDA,  Q,  CLAR,  E0TA,  SW,  N0EIGJ,  NS, 
ICALCS,  IPRINS,  IPRINA) 

Definition  of  call  statement  variables: 


( ) ■ program  dimensions 


Variable  ( ) 

Data  TVpe 

Description 

AFW(10,3) 

INPUT 

Rigid  load  distributions  on  NS  number  of 
chordwise  strips.  Stored  by  columns: 

Col  1,  a load;  col  2,  6 FLAP  loacl‘> 
col  3,  inertia  load  for  Nz  * 1. 

AFWCP(10,3) 

INPUT 

Chordwise  centers  of  pressure  for  rigid 
loads  provided  in  AFW.  These  X/C  values 
are  stored  by  columns. 

AFWF(10-3) 

OUTPUT 

Total  flexible  load  distributions,  by 
column. 

AFWDF(10,3) 

OUTPUT 

Delta  flexible  load  distributions,  by 
column. 

F(3,4) 

OUTPUT 

Simulation  of  load  distributions:  Col  1, 
I rigid  loads;  col  2,  Z total  flexible 
loads;  coi  3,  Z delta  flexible  loads; 
col  4,  flexible-to-rigid  ratios. 

YEIGJ(20) 

INPUT 

Elastic  axis  y-coordinates  for  wing  El 
and  GJ.  NOTE:  YEIGJ  (1)  must  equal 
zero,  (in.). 

El (20) 

INPUT 

2 

El,  wing  bending  stiffness  (lb  in.). 

GJ(20) 

INPUT 

2 

GJ,  wing  torsional  stiffness  (lb  in.). 
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Variable  ( ) 

Data  Type 

Descriptions 

CR 

INPUT 

Chord  of  wing  at  fuselage  <£  (in.). 

602 

INPUT 

Wing  semispan,  b/2  (in.). 

BIBS 

INPUT 

Butt  line  of  body  side  (in.). 

ANGLE 

INPUT 

Angle  of  wing  leading  edge  sweep  (deg) 

FSLERT 

INPUT 

Fuselage  station  of  wing  leading  edge 
at  fuselage  £ (in.). 

ANGEA 

INPUT 

Angle  of  wing  elastic  axis  sweep  (deg) 

FSEA0 

INPUT 

Fuselage  station  of  elastic  axis  at  <£ 
(in.). 

WLAMDA 

INPUT 

Wing  taper  ratio, 

Q 

INPUT 

Dynamic  pressure , q (psf ) . 

CLAR 

INPUT 

CLft  /RAD  for  the  wing. 

E0TA 

INPUT 

Exposed/total  load  ratio,  for 
distributions. 

SW 

INPUT 

Wing  area  (ft^). 

N0EIGJ 

INPUT 

Number  of  values  in  the  YEIGJ,  El  and 
GJ  lists,  20  maximum. 

NS 

INPUT 

Numbers  of  load  strips  on  the  exposed 
wing,  10  maximum. 

ICALCS 

INPUT 

Controls  calculation  of  SIC  matrix: 

0 - Do  not  calculate  SIC 

1 - Do  calculate  SIC 
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Variable  ( ) 


Data  Type 


Descriptions 


IPRINS 


IPRINA 


INPUT  Controls  printing  of  SIC  matrix  and 
other  diagnostic  data  from  subroutine 
FLXSIC: 

0 - Do  not  print  SIC 

1 - Do  print  SIC 

NOTE:  IPRINS  has  no  effect  if  ICALCS  - 0 

INPUT  Controls  printing  of  diagnostic  data 
from  subroutine  WFLEX: 

0 - Do  not  print  diagnostic  data 

1 - Do  print  diagnostic  daca 


Restrictions 


The  number  of  load  strips,  NS,  may  not  exceed  10,  3nd  the  number  of  El 
and  GJ  values,  N0EIGJ,  may  not  exceed  20. 


Error  Returns 


If  NS  is  greater  than  10,  or  N0EIGJ  is  greater  than  20,  an  error 
message  will  be  printed  and  execution  will  be  stopped. 


Method 


The  method  is  described  under  the  heading  "Wing  Aeroelastic  Distributions" 
in  Section  II.  Note  that  when  using  subroutine  GLSQ  to  obtain  the  solution, 
the  D-matrix  must  be  reloaded  into  A for  each  load  effect,  because  GLSQ 
destroys  A in  the  solution  process.  Also,  note  that  region  A must  be  cleared 
prior  to  loading  D for  at  least  two  more  rows  than  the  number  of  rows  in  D, 
in  order  for  the  GLSQ  solution  to  work. 
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Operational  Notes 


Rigid  load  distributions  that  enter  into  WFLEX  are  for  loads  on  the 
exposed  wing  panel  only  and  can  be  handled  several  ways.  But  regardless  of 
the  manner  in  which  it  is  used,  the  a load,  per  radian,  on  the  exposed  wing 
panel  due  to  a unit  angle  of  attack  is  computed  in  the  progran  as: 

NS 

LOAD  - CLAR*EjjrA*Q*^*  £ AW(I,1) 

1-1 

The  variable  E0TA  in  application  to  this  flexible  airloads  program 
is  equated  to  1.0. 


SUBROUTINE  FLXSIC 

Using  wing  geometry  and  wing  El  and  GJ  distributions  as  input  data, 
subroutine  FLXSIC  computes  wing  structural  influence  coefficient  for  use 
in  the  static  aeroelastic  solution  that  is  calculated  in  subroutine  WFLEX. 
The  logical  flow  chart  and  program  listing  are  shown  in  Appendix  A. 


Usage 


All  input-output  data,  except  printed  diagnostic  results,  are  passed 
through  the  call  statement. 

CALL  FLXSIC  (SICBAR,  YEIGJ,  El,  GJ,  CR,  B02,  BLBS,  ANGLE,  FSLERT, 
ANGEA,  FSEA0,  WLAMDA,  X0CFWD,  X0CAFT,  N0EIGJ,  NS,  IPRINS) 


Definition  of  call  statement  variables: 
( ) * Program  dimensions 


Variable  ( ) Data  Type 


Description 


SICBAR(10,20) 


OUTPUT  Matrix  of  strip  slopes  due  to  loads  at 
the  SIC  points  (radians/lb). 


YEIGJ(20)  INPUT  Elastic  axis  y-coordinates  for  wing  El 

and  GJ.  NOTE:  YEIGJ (1)  must  equal 
zero,  (in.). 
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Variable  ( ) 

Data  Type 

Description 

El (20) 

INPUT 

2 

El,  wing  bending  stiffness  (lb-in.) • 

GJ(20) 

INPUT 

2 

GJ,  wing  torsional  stiffness  (lb-in.). 

CR 

INPUT 

Chord  of  wing  at  fuselage  (£  (in.). 

B02 

INPUT 

Wing  semispan,  b/2  (in.). 

BLBS 

INPUT 

Butt  line  of  body  side  (in.). 

ANGLE 

INPUT 

Angle  of  wing  leading  edge  sweep  (deg). 

FSLERT 

INPUT 

Fuselage  station  of  wing  leading  edge 
at  fuselage  <£  (in.). 

ANGEA 

INPUT 

Angle  of  wing  elastic  axis  sweep  (deg). 

FSEA0 

INPUT 

Fuselage  station  of  elastic  axis  at  G 
(in.). 

WLAMDA 

INPUT 

Wing  taper  ratio,  ^ * “P^d 

X0CFWD 

INPUT 

Strip  X/C  forward,  used  as  0.1S. 

X0CAFT 

INPUT 

Strip  X/C  aft,  used  as  0.65. 

N0EIGJ 

INPUT 

^ber  of  values  in  the  YEIGJ,  El  and 
GJ  lists,  20  maximum. 

NS 

INPUT 

Number  of  load  strips  on  the  exposed 
wing,  10  maximum. 

IPRINS 

INPUT 

Controls  printing  of  SIC  matrix  and 
other  diagnostic  data: 

0 - Do  not  print  SIC  and  diagnostic 

1 - Do  print  SIC  and  diagnostic 
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Restrictions 


The  number  of  load  strip,  NS,  may  not  exceed  10,  the  lumber  of  El  and 
GJ  values  may  not  exceed  20,  and  YEIGJ(l)  mist  equal  0. 


Error  Returns 

In  the  logic  of  this  routine,  it  is  possible  for  y in  the  swept  elastic 
axis  system  to  become  negative  for  a highly  swept  wing.  If  it  does,  this 
subroutine  will  print  out  an  error  message  and  execution  will  be  stopped. 

Method 

The  method  is  described  under  the  heading  "Wing  Structural  Influence 
Coefficients,"  in  Section  II.  Coding,  with  the  exception  of  the  technique 
used  in  integration  of  equation  44  is  self-explanatory.  Integration  of 
equation  44  is  accomplished  trapezoidally  after  the  input  El  and  GJ  data 
are  first  interpolated,  using  the  C0DIM2  subroutine,  to  obtain  20  intervals 
between  the  £ and  the  point  at  which  the  deflection  is  being  computed.  For 
example: 

fy  a 

^0  i«2 

where 

AX  « y/20 

Note  that  if  the  value  of  y ever  exceeds  the  value  of  n , the  integration 
upper  limit  is  changed  to  "n  by  simply  zeroing  out  the  value  of  the  func- 
tions of  X.  in  both  the  bending  and  torsion  integrals. 
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Recommended  Operational  Data  Setup 


In  the  formulation  of  the  aeroelastic  solution  that  uses  the  SIC  matrix 
from  this  subroutine,  the  wing  elastic  axis  (EA)  has  been  carried  inboard  of 
the  side  of  the  body  to  the  body  centerline.  In  general,  the  wing  carry- 
through  structure  does  not  extend  along  the  wing  EA,  but  is  directed  across 
the  airplane.  Since  bending  across  the  airplane  does  not  affect  wing  angle  of 
attack,  and  the  wing  is  effectively  fixed  in  torsion  at  the  sire  of  the  body, 
stiffness  El  and  GJ  inboard  of  the  body  side  should  be  arbitrarily  greatly 
increased. 

El  and  GJ  recommended  data  setup: 

BSs 

XBODY  SIDE  “ cos  A c 

eA 


YEIGJ(l) 

■ 0. 

EI(1) 

- 10 

eibs 

Gj(l) 

- 10 

"bs 

(2) 

- 0.25 

^BS 

(2) 

- 10 

eibs 

(2) 

» 10 

“bs 

(3) 

* 0.50 

^BS 

(3) 

- 10 

eibs 

(3) 

- 10 

^BS 

(4) 

" ^BS 

(4) 

*eibs 

(4) 

* GJ_„ 
BS 

(5) 

« 1.10 

7BS 

(5) 

* See  Note 

(5) 

« See  Note 

NOTE:  Data  for  point  5 and  on  should  be  taken  from  the  El  and  GJ  curves. 

Sufficient  data  points  should  be  input  (20  maximun)  to  give  good  defi- 
nition of  the  curves  when  fitted  with  C0DIM2. 

Subroutines  USPANF  through  RERDAT  set  up  El  and  GJ  values  in  accordance 
with  the  preceding  note  and  interpolate  El  and  GJ  input  data  BF(47)  through 
BF(106)  to  obtain  15  additional  points  outboard  of  YEIGJ(5).  Resulting  El 
and  GJ  data  are  then  located  in  BF(114)  through  BF(173) , as  shown  in  Table  4, 
for  transfer  to  this  subroutine. 


SUBROUTINE  GLSQ 

This  FORTRAN  subprogram  gives  the  least  squares  solution  to  a system  of 
overdetermined  linear  equations  Bx  ■ C,  where  B is  an  N x M matrix,  with 
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N » M and  C a colimn  vector  of  dimension  N.  The  logical  flow  chart  and  the 
program  listing  are  shown  in  Appendix  A. 


Calling  Sequence 

CALL  GLSQ  (A,  X,  IL,  N,  M,  ALPHA,  El,  E2) 
Definition  of  call  statement  variables: 


( ) - program  dimensions 


Variable  ( ) 

Data  Type 

Description 

A(25,26) 

INPUT 

The  augmented  matrix  B,C  of  at  least 
dimension  (M+2)  x (M+l)  or  N x (M+l) , 
whichever  is  greater. 

X(26) 

OUTPUT 

Die  vector  where  the  solution  is  stored 
and  must  be  of  at  least  dimension  M+  1. 

IL (26) 

INPUT 

A temporary  vector  of  at  least  dimen- 
sion M+l. 

N 

INPUT 

Number  of  rows  of  B. 

M 

INPUT 

Number  of  colunns  of  B. 

ALPHA 

OUTPUT 

The  square  root  of  the  sum  of  the 
squares  of  the  residuals. 

El  and  E2 

INPUT 

Two  nonnegative  numbers  which  are 
small  compared  to  the  size  of  the  num- 
bers in  the  input  matrix. 

Restrictions 

This  subprogram  has  a dimension  statement  A(25,  26),  X(26),  IL(26). 

The  calling  program  must  be  dimensioned  exactly  the  same.  The  input  data,  A, 
are  destroyed  during  the  computation. 
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Method 


The  method  is  discussed  under  "Wing  Aeroelastic  Distributions,"  in 
Section  II.  In  the  triangularization  process,  if  any  number  to  be  annihilated 
is  already  less  in  magnitude  than  El,  it  is  considered  to  be  zero  and  the 
computations  to  annihilate  the  element  are  omitted. 

Let  bj  denote  the  column  of  3.  If  the  routine  finds  that  bk  is  a 
linear  combination  of  bj,...,  b^.^,  then  it  sets  the  k*h  component  of  x to 
zero,  The  routine  considers  that  such  a linear  combination  holds  whenever 
lumbers  X},....,  Xk-i  can  be  found  so  that  the  length  of  the  vector 


V 


(Xjbj  +....♦  Vk_j  bk_j) 


is  less  than  E2. 


Operational  Notes 

When  using  GLSQ  to  solve  an  exactly  determinate  set  of  linear  equations, 
the  dimension  of  the  augmented  matrix  must  be  at  least  A(M+2,  M+l) , and  the 
region  must  be  cleared  prior  to  loading  the  augmented  matrix. 


SUBROUTINE  MATRIT 

Subroutine  MATRIT  prints  matrices  by  row  or  column  in  format  IP8E13.5 
and  with  a heading  statement  that  is  transmitted  in  the  calling  sequence. 
The  logical  flow  chart  and  the  program  listing  are  shown  in  Appendix  A. 


Usage 


CALL  MATRIT  (AMN,  NR,  NC,  NRMAX,  MTYPE,  IPRIN,  HEAD) 
Definition  of  call  statement  variables: 

( ) = Program  Dimensions 
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Variable  ( ) 

Data  Type 

Description 

AMN(l) 

INPUT 

Matrix  to  be  printed. 

NR 

INPUT 

Number  of  rows  in  AMN. 

NC 

INPUT 

Number  of  colurai  in  AMN. 

NRMAX 

INPUT 

Maximun  number  of  rows  in  AMN, 
dimension  in  calling  program. 

NfTYPE 

INPUT 

Type  of  matrix: 

1 - real 

-1  - complex  (do  not  use) 

IPRIN 

INPUT 

Print  by  rows  or  columns: 

1 - print  by  rows 

2 - print  by  columns 

HEAD(IOO) 

INPUT 

Format  of  matrix  heading;  for 
example , 

24H  (IHI  20X,  11H4ASS,  MATRIX). 

SUBROUTINE  BNLDSF 

Subroutine  BNLDSF  is  used  to  determine  the  limit  airloads  and  CP's  an 
the  airplane  components,  and  up  inertia  factors  for  the  flight  conditions 
specified  by  the  main  program.  There  are  five  types  of  flight  conditions  for 
which  the  airloads  can  be  determined;  namely,  (1)  the  balanced  maneuver  con- 
dition, (2)  the  pitching  acceleration  condition,  (3)  the  vertical  gust  condi- 
tion, (4)  the  lateral  gust  condition,  and  (5)  the  yawing  acceleration 
condition.  The  methods  used  in  deriving  the  component  loads  for  these 
conditions  are  presented  in  Section  II  and  include  the  aeroelastic  effects  of 
wing  flexibility. 

The  subroutine  BNLDSF  logical  flow  chart  and  the  program  listing  are 
shown  in  Appendix  A. 
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SUBROUTINE  SPABMF 


Subroutine  SPABMF  is  used  to  determine  the  net  limit  airload  shear  and 
moments  along  the  load  reference  line  on  the  wing,  horizontal  tail,  and  verti- 
cal tail  surfaces.  The  methods  employed  are  described  in  Section  II.  The 
logical  flow  chart  and  the  program  listing  are  shown  in  Appendix  A. 


OPERATING  CORE,  TC0M 

The  program  operating  core,  TO0M,  is  dimensioned  to  4,400  locations  and 
is  blocked  into  data  regions  as  shown  in  Table  10.  A description  and  the 
location  of  the  data  within  each  data  block  can  be  found  in  the  tables 
referenced  by  Table  10. 


OUTPUT  DATA  DESCRIPTION 
GENERAL 

Primary  output  of  the  program  is  a deck  of  cards  containing  airplane 
component  limit  airloads,  CP  locations,  and  wing  and  empennage  limit  airload 
shears  and  moments  at  the  spanwise  stations  along  the  load  reference  line 
selected  for  weight  analysis.  Punched  card  output  is  in  a format  for  use  as 
an  optional  airloads  input  to  the  SWEEP  program.  Printout  of  other  data  is 
also  included  for  visual  inspection  of  the  results  of  final  and  intermediate 
calculations. 


PUNCHED  CARD  OUTPUT 

The  output  punched  cards  list  limit  airloads  data  in  an  E-foimat  to  be 
read  by  subroutine  DECRD  in  SWEEP.  A sample  printout  of  data  on  the  punched 
cards  is  shown  in  Table  17.  Each  card  contains  a DECRD  array  index  number, 
followed  by  five  decimal  data  items  and  the  card  identification  number. 

A "yes"  in  control  factor  ND(4Q),  per  Table  1,  will  provide  the  punched  card 
output . 

Output  data  array  is  described  in  Table  15.  The  index  number  punched  on 
each  card  corresponds  to  the  location  within  the  array  for  the  first  item  of 
decimal  data  on  the  card. 
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Condition  identification  nunber  and  component  loads  identification 
nunbers  are  assigned  by  the  program  and  are  included  in  the  punched  output 
deck  for  each  condition.  Component  loads  identification  nunbers  also  appear 
in  the  printed  output. 


Condition  ID 


The  condition  identification  nunber  is 
card  1 and  is  defined  as  follows: 


punched  in  columns  13  to  24  of 


f 


Case  number 

(Refer  to  table  2) 


Condition  array  ID 
(Refer  to  table  18) 


W V 


Wing  ID 

1 - Fixed 

2 - VS  - fwd 

3 - VS  - aft 


I 

Vertical  tall  ID 

5 - Single  tali 

6 - Dual  tai  1 

7 - T-tall 


Fuselage  Loads  ID 


The  fuselage  loads  identification  number  is  punched  in  columns  25  to  36 
of  card  1 and  appears  on  the  components  load  and  CP  printed  page.  It  is 
defined  as  follows: 

CA  CN  0 


Wingloads  ID 


The  wingloads  identification  nunber  is  punched  in  columns  25  to  36  of 
card  5 and  appears  on  the  wing  spanwise  loads  printed  page.  It  is  defined  as 
follows: 


CA  CN  W 
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Horizontal  Tail  Loads  ID 


The  horizontal  tailloads  identification  nunber  is  punched  in  coluims  61 
to  72  of  card  15  and  appears  on  the  horizontal  tail  spanwise  loads  printed 
page.  It  is  defined  as  follows: 

CA  CN  4 


Vertical  Tail  Loads  ID 

The  vertical  tailloads  identification  nunber  is  punched  in  columns  37 
to  48  of  card  26  and  appears  on  the  vertical  tail  spanwise  loads  printed  page. 
It  is  defined  as  follows: 

CA  CN  V 


Punched  Card  ID 

The  punched  output  identification  and  card  sequence  nunber  is  punched  in 
coluims  73  to  80  and  is  defined  as  follows: 


1 CA  CN 


Signifies  sweep  1 data 


SEQ 


Three-column  sequence  No. 


PRINTED  OinWT  DATA 

The  program  will  always  print  the  data  in  Tables  19  through  28.  With  a 
"yes"  in  ND(41),  per  Table  1,  a listing  of  the  data  on  the  punched  output 
cards  will  be  printed  as  shown  in  Table  17.  With  a "yes"  in  ND(44) , per 
Table  1,  FLXSIC  matrices  will  be  printed  as  shown  in  Tables  29  through  35. 
With  a "yes"  in  ND  (45),  per  Table  1,  WFLEX  matrices  will  be  printed  as  shown 
in  Tables  36  through  38. 
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Data  of  Primary  Interest 


Data  of  primary  interest  to  the  user  consist  of  (1)  that  which  is 
descriptive  of  the  data  that  are  to  be  used  as  optional  input  to  SWEEP,  and 
(2)  input  data  for  the  main  program,  BFCNTL. 

Table  17  shows  a sample  printout  of  the  data  on  the  output  punched  cards. 
Data  items  can  be  identified  using  Table  15.  Table  25  presents  a summary 
of  the  component  limit  airloads,  CP's,  and  the  airplane  inertia  factors. 

Tables  26  through  28  show  sample  printouts  of  limit  airload  shears  and 
moments  along  the  load  reference  line  for  the  wing,  horizontal  tail,  and 
vertical  tail,  respectively.  Tables  19  and  20  show  sample  printouts  of  the 
program  input  BC  and  BF  data  sets;  the  data  items  can  be  identified  using 
Tables  3 and  4. 


Intermediate  Step  and  Diagnostic  Data 


Tables  of  printed  data  presented  are  divided  into  three  groups.  The 
first  group  consists  of  intermediate-step  results  produced  by  subroutine 
USPANF  as  spanwise  unit  airload  distributions  for  the  wing,  and  horizontal  and 
vertical  tails  are  developed.  The  second  group  consists  of  optional  diag- 
nostic matrix  data  produced  by  subroutine  FLXSIC.  The  third  group  consists  of 
optional  diagnostic  matrix  data  produced  by  subroutine  WFLEX. 


Intermediate-Step  Data 

Tables  21,  23,  and  24  are  produced  by  subroutine  USPANF.  Table  22  is 
produced  by  subroutine  WFLEX,  which  is  subordinate  to  subroutine  USPANF. 

These  printed  data  are  not  optional  in  the  current  configuration  of  the  flexi- 
ble airloads  stand-alone  program.  Table  21  shows  sample  printouts  of  the 
rigid  loading  data  extracted  from  the  data  bank  aerodynamic  data  for  the 
condition  mach  number  and  the  subject  air  vehicle  configuration.  Table  22 
shows  sample  printouts  of  r’gid  loading  and  CP  data  inputed  to  subroutine 
WFLEX  by  subroutine  USPANF,  and  aeroelastic  loading  data  returned  to  sub- 
routine USPANF  by  subroutine  WFLEX.  Table  23  shows  sample  printouts  of  final 
unit  spanwise  loading  distributions  3 or  the  wing,  and  horizontal  and  vertical 
tails  developed  by  subroutine  USPANF.  Table  24  shows  sample  prinumts  of 
spanwise  variation  of  unit  shear,  bending  moment,  and  torque  for  three  sur- 
faces as  developed  by  subroutine  USPANF  and  stored  in  data  region  BU. 
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FLXSIC  Diagnostic  Data 

Tables  29  through  35  present  optional  diagnostic  data  produced  by 
subroutine  FLXSIC.  If  control  factor  ND(44)  is  assigned  a value  of  1,  these 
tables  will  be  produced.  A value  of  0 will  eliminate  them.  Table  29  shows 
a sample  printout  of  the  data  inputed  to  subroutine  FLXSIC  by  subrou- 
tine WFLEX.  Table  30  shows  a sample  printout  of  the  elastic  axis  geometry 
developed  by  subroutine  FLXSIC.  Tables  31  through  35  show  sample  printouts  of 
the  various  structural  influence  coefficient  (SIC)  matrices  developed  by 
subroutine  FLXSIC.  Subroutine  FLXSIC  is  recalled  only  when  there  is  a change 
in  wing  geometry  or  stiffness  properties. 


WFLEX  Diagnostic  Data 

Tables  3 b through  38  present  optional  diagnostic  data  produced  by  sub- 
routine WFLEX.  If  control  factor  ND(45)  is  assigned  a values  of  1,  these 
data  will  be  produced.  A value  of  0 will  eliminate  them.  Tc'  le  36  shows  a 
sample  printout  of  the  pertinent  aerodynamic  data  used  by  subroutine  WFLEX. 
Table  37  shows  a sample  printout  of  the  aeroelastic  D-matrix  developed  by  the 
subroutine.  Table  38  shows  sample  printouts  of  the  rigid  and  the  resultant 
flexible  loads  matrices.  Subroutine  WFLEX  is  called  whenever  mach  number, 
dynamic  pressure,  wing  geometry,  wing  stiffness  properties,  or  wing  weight 
distribution  change. 
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TABLE  1.  ND  CONTROL  FACTORS  CARD 


Control  Factor 


SWEEP  A/V  class 
Fighter 
Attack 
Bomber  I 
Bomber  II 
Cargo-assault 
Cargo- transport 


Wing  Type 
Fixed  wing 
Variable  sweep 


Vertical  tail  type 
Single  tail 
Dual  tail 
T-tail 


Loads  requirements 
Select  conponents 
Fuselage  No 

Yes 

Wing  No 

Yes 

Horizontal  tail  No 
Yes 

Vertical  tail  No 
Yes 

Do  all  conponents 


ND  ND  ND  ND  ND 

13  14  15  23  24 


ND  ND  ND 

25  26  27 


22,24 
26,28 
5 30 
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TABLE  1.  ND  CONTROL  FACTORS  CARD  (CONI’) 


Control  Factor 


Basic  condition  types* 
Pos  bal  flight 

Neg  bal  flight 

Maneuvering  flap 

1 g trim  flap 

Pos  vert  gust 

Neg  vert  gust 

Lateral  gust 

Pitch  acceler 

Yawing  acceler 


runch  output  factors 
Punched  output 

Print  check 


Card  ND  ND  ND  ND  ND  ND  ND  ND  ND  ND  ND 

Col  28  29  30  31  32  33  34  35  36  40  41 


No 
Yes 
No 
Yes 
No 
Yes 
No  37,38 

Yes 38 

No  39,40 

Yes 40 

No  41,42 

Yes 42 

No  43,44 

Yes 44 

No  45, 
Yes 


No  55,56 

Yes 56 

No  57,58 
Yes  58 
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TABLE  1.  ND  CONTROL  FACTORS  CARD  (GONCL) 


Control  Factor 

Card 

Col 

ND 

42 

ND 

43 

ND 

44 

ND 

45 

ND 

46 

Input  control  factors  - 
No.  of  values 
YEIGJ,  EL,  GJ 
(20  Maximum) 

Wing  weight  strips 

■ 

20 

61,62 

10 

Diagnostic  printing 

FLXSIC  matrices 

No 

64 

0 

Yes 

64 

1 

MFLEX  matrices 

No 

66 

0 

Yes 

66 

1 

Case  No. 

67,68 

** 

*Basic  conditions  not  specifically  defined  by  the  chosen  case  nunber  will 
not  be  produced,  even  if  the  type  factor  is  entered  as  yes.  Type  factors 
are  used  for  the  selection  or  rejection  of  basic  conditions  described  by 
the  chosen  case  number. 

**Enter  case  number  selected  from  Table  2 (the  number  must  be  right- 
adjusted  in  the  field). 
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TABLE  2.  CASE  NUMBER  CONDITIONS 


Case  No.  Condition  No. 


108  WV 

110  WV 

111  WV 

112  WV 

114  WV 

115  WV 


208 
210 

211  WV 

212  WV 

214  WV 

215  WV 


315  WV 


401  WV 

402  WV 

403  WV 
405  WV 

409  WV 
413  WV 


509  WV 
513  WV 


Condition 

Type 

Mach  No. 

+ Bal  man. 
Pitch  accel 
Yaw  accel 

+ Bal  man. 
Pitch  accel 
Yaw  accel 

Mli,i 

♦ Bal  man. 
Pitch  accel 
Yaw  accel 

+ Bal  man. 
Pitch  accel 
Yaw  accel 

mL1,2 

+ Bal  man. 
Pitch  accel 
Yaw  accel 

+ Bal  man. 
Pitch  accel 
Yaw  accel 

^1,3 

+ Vert  gust 

- Vert  gust 
Lat  gust 

- Bal  man. 

Vi 

- Bal  man. 

- Bal  man. 

+ Vert  gust 
- Vert  gust 
Lat  gust 

Mh1,2  or 

20,000  ft* 

- Bal  man. 

- Bal  mar. 

MH1 , 2 

Altitude 
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TABLE  2.  CASE  NUMBER  CONDITIONS  (CONCL) 


Condition 

Case  No. 

Condition  No. 

Type 

Mach  No. 

Altitude 

601  WV 

+ Vert  gust 

*k,3  or 

£ 

602  WV 

- Vert  gust 

20,000  ft* 

h,  . or  20,000  ft* 

0 

603  WV 

Lat  gust 

1 9 3 

613  WV 

- Bal  man. 

Mm, 3 

hl*3 

1 

+ Bal  man. 

- Bal  man. 

hl  1 

Pitch  man. 

707  WV 

Yaw  man. 

708  WV 

+ Bal  man. 

*7 

709  WV 

- Bal  man. 

0.90 

7 

710  WV 

Pitch  man. 

hl,2 

711  WV 

Yaw  man. 

712  WV 

+ Bal  man. 

713  WV 

- Bal  man. 

U 

714  WV 

Pitch  man. 

hl,3 

715  WV 

Yaw  man. 

8 

816  WV 

Flap  man. 

1.5  Vso 

SL 

9 

917  WV 

1 g flap  trim 

1.2  VSL 

SL 

1018  WV 

+ Vert  gust 

1019  WV 

- Vert  gust 

^2  1 

10 

1020  WV 

+ Bal  man. 

Vi 

L,l 

1021  WV 

- Bal  man. 

1022  WV 

+ Bal  man. 

Vi 

1023  WV 

- Bal  man. 

h2,2 

11 

1118  WV 

1119  WV 

+ Vert  gust 
- Vert  gust 

Mr2,2  or 
Mh2,  20,000  ft* 

h 0 or  20,000  ft* 

2,4. 

1122  WV 

+ Bal  man. 

^2,2 

Vi 

1123  WV 

- Bal  man. 

h2,2 

*Gust  conditions  maxinum  altitude  is  20,000  feet. 


W ■ 1 « Fixed  wing  V ■ 5 * Single  vertical  tail 

2 ■ Variable  sweep  (fwd)  6 ■ Dual  vertical  tail 

3 - Variable  sweep  (aft)  7 * T - tail 
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TABLE  3.  BC  INPUT  DATA  SET 
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28-30  SPEEDS  MH  FOR  VARIABLE  SHEEP  HI  NG  IN  FOPWARD  POSITION.  MACH  NUMBER 

31  MINIMUM  SPEED  FLAPS  UP  IVSOI  AT  MDH.  KNOTS 

32  MINIMUM  SPEED  FLAPS  DOUN  ( VSL ) AT  LDH.  KNOTS 

33  XO— DI STANCE  FROM  X REFERENCE  POINT  TO  BODY  NOSE.  INCHES 

34  LN-LENGTH  OF  BODY  NOSE.  INCHES 


TABLE  3.  BC  INPUT  DATA  SET  (OONT) 
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TABLE  5.  DT  DATA  SET 


Location 

Description 

DT(1)  is  equivalent  to  TQZM  (100) 

1 - 4 

(1)  through  (4),  span  station,  fraction  of  semispan 

5 - 9 

(1)  through  (4),  taper  ratios 

9 - 20 

M(l)  through  M(12),  mach  numbers 

21  - 32 

B/K  (1)  through  B/K  (12),  compressible  section  lift- 
curve  slope  parameters 

33  - 44 

X/C  (1)  through  X/C  (12) , wing  section  center  of 
pressure  locations,  wfraction  of  chord 

45  - 56 

X/Gp  (1)  through  X/Gp  (12),  tail  section  center  of 
pressure  locations,  fraction  of  chord 
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TABLE  6.  DB  DATA  SET 


Location  Description 

DB(1)  is  equivalent  to  TG0M  (156) 


1-6 

Ag(l)  through  Ag(6),  compressible  sweep  angle  (deg) 

7 - 13 

BA/K(1)  through  BA/K(6),  aspect  ratio  parameter 

14  - 55 

Table  of  C | C/ClC^v  values  versus  Ag  and  BA/K  for  n 1) , 

\(1)  (42  values  per  table) 

56  - 97 

Same  as  14  - 55,  except  for  n(l),  \(2) 

78  - 139 

Same  as  14  - 55,  except  for  n(l)»  *-(3) 

140  - 181 

Same  as  14  - 55,  except  for  n(l),  X(4) 

182  - 349 

Same  as  14  - 181,  except  n * n(2) 

350  - 517 

Same  as  14  - 181,  except  n ■ -n  (3) 

518  - 685 

Same  as  14  - 181,  except  n ■ n(4) 

686  - 727 

Table  of  BCl^/K  values  versus  Ag  and  BA/K  for  \(1) 

(42  values  per  table) 

728  - 769 

Same  as  686  - 727,  except  \«  \(2) 

770  - 811 

Same  as  686  - 727,  except  \ - \(3) 

812  - 853 

Same  as  686  - 727,  except  K-  X(4) 

NOTES:  Table  data  order  ■ values  for  BA/K(1),  Ag(l)  through  Ag(6) 

BA/K(2),  Ag(l)  through  Ag(6),  etc 

Refer  to  Table  2 in  Volume  III,  Section  II,  for  numerical  values  of 

data. 
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TABLE  7.  DF  DATA  SET 


Location 

Description 

DF(1)  is  equivalent  to  TO0M  (1009) 

1 - 5 

Cp/%(1)  through  Cp/Cw(5),  flap  chord  to  wing  chord  ratios 

6 - 10 

X/C  (1)  through  X/C  (5),  centers  of  pressure  of  section 
lift  due  to  flap  deflection,  fraction  of  wing  chord 

11  - 15 

Kq:  (1)  through  K^f  (5),  flap  lift  effectiveness 
parameters 

16  - 25 

bp/btf  (1)  through  bp/bw  (10) , flap  span  to  wingspan 
ratios 

26  - 36 

t\  (1)  through  "n (11)  wingspan  stations,  fraction  of 
wing  semispan 

37  - 146 

Table  of  Cl  C/CLtf  C^y  values  versus  n , bp/b^ 

NOTE:  Table  data  order  ■ values  for  bp/b^l),  n (1)  through  n(ll) 

bp/b^),  n(l)  through  n(ll),  etc 


Refer  to  Table  3 in  Volume  III,  Section  II,  for  numberical  value  of 

data. 
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TABLE  8.  DP  DATA  SET 


Location 

Description 

DP(1)  is  equivalent  to  TO0M  (1155) 

1 - 5 

BA(1)  through  M(5),  aspect  ratio  parameter 

6 - 14 

fim(l)  through  Bm(9) , sweep  parameter 

15  - 59 

Table  of  2r/Vbb  values  versus  Bm  and  BA  for  n (1), 
\ (1)  (45  values  per  table) 

60  - 104 

Same  as  15  - 59,  except  for  n(l)»  X.  (2) 

105  - 149 

Same  as  15  - 59,  except  for  n(l)»  X(3) 

150  - 194 

Same  as  15  - 59,  except  for  n(l)»  K(4) 

195  - 374 

Same  as  15  - 194,  except  ■ n(2) 

375  - 554 

Same  as  15  - 194,  except  tj  ■ n(3) 

555  - 734 

Same  as  15  - 194,  except  n-  n(4) 

NOTE:  Table  data  order  • Values  for  BA(1),  Bm(l)  through  Bm(9),  etc 

Refer  to  Table  4 in  Volume  III,  Section  II,  for  numerical  values  of 

data. 
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DATA  LOCATIONS  87-128  ARE  FOR  A VARIABLE  SWEEP  WING  IN  THE  FORWARD  POSITION. 

87  WING  LEADING  EDGE  SWEEP  ANGLE.  OEGREES 

88  WING  ELASTIC  AXIS  SWEEP  ANGLE.  OEGREES 

89  WING  LOADS  REFERENCE  AXIS  SWEEP  ANGLE.  DEGREES 
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TABLE  10.  TO0M  BLOCKED  CORE  LOCATIONS 
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TABLE  11.  CONDITION  DATA  REGION  BB*** 


Location 


Description 


BB(1)  is  equivalent  to  TC0M(2953) 


W,  airplane  weight  (lb) 

XCG,  airplane  center  of  gravity  location  (FS)  (in.) 

Iy,  airplane  pitching  moment  of  inertia  (slug  ft?) 
lit  airplane  yawing  moment  of  inertia  (slug  ft2) 

6p,  flap  deflection  (deg) 

N2,  airplane  limit  normal  load  factor 

Q,  airplane  pitching  acceleration  (rad/sec2) 

R,  airplane  yawing  acceleration  (rad/sec2) 
hp,  altitude  (ft) 

M,  mach  number 

Previous  condition  altitude  (ft) 

Wing  ID  number* 

Horizontal  ID  number  "4.0 
Vertical  tail  ID  number** 

Condition  number 

Previous  condition  mach  number 

q,  dynamic  pressure  (psf) 

cLoWR»  wing  rigid  lift  curve  slope  per  radian 

Vg,  equivalent  airspeed  (knots) 

Balance  , airplane  normal  load  factor  for  balanced  part 


*Fixed  wing  ID  ■ 1.0 
Variable  wing  forward  ID  ■ 2.0 
Variable  wing  aft  ID  ■ 3.0 

**Single  vertical  tail  ID  - 5.0 
Dual  vertical  tail  ID  - 6.0 
T-type  vertical  tail  ID  - 7.0 

***This  condition  data  region  is  set  up  for  each  condition 
by  the  main  program  BFCNTL. 
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TABLE  12.  PARAMETER  DATA  REGION  BS  FOR  SUBROUTINES 
USPANF  AND  BNLDSF 


Location 

Description 

USPANF 

BS  (I ) is  equivalent  to  TQ0M(2973) 

1 

B/K,  subsonic  lift-curve  slope  parameter 

2 

M,  mach  number 

| 3 

(X/C)cp,  section  center  of  pressure,  fraction  of  chord 

4 

Ygg»  wing-body  interface  station  (BL),  or  horizontal 

5 

Ale,  surface  leading  edge  sweep  angle  (deg) 

6 

A£A,  surface  elastic  axis  sweep  angle  (deg) 

7 

A pa,  surface  load  reference  line  sweep  angle  (deg) 

8 

Xle»  fuselage  sta.  of  root  chord  leading  edge  (in.) 

9 

Xea,  fuselage  sta  of  elastic  axis  at  root  chord  (in.) 

10 

Cr,  root  chord  of  theoretical  surface  (in.) 

11 

\s,  taper  ratio  of  theoretical  surface 

12 

Ag,  aspect  ratio  of  theoretical  surface 

13 

Ss,  area  of  theoretical  surface  (ft2) 

14 

bg/2  or  by,  span  of  surface  panel  (in.) 

15 

nBS  ornBV»  sVm  station  of  surface-body  interface 
station,  fraction  of  panel  span 

16 

(X/C)ra,  location  of  load  reference  line  in  fraction  of 
chord 

17 

Ypo/bw/2»  wingspan  station  of  outboard  end  of  flap, 
fraction  of  semispan 

18 

Ypi/tty/2,  wingspan  station  of  inboard  end  of  flap, 
fraction  of  semispan 

19 

Ab»  subsonic  compressible  sweep  parameter,  or  BA, 
supersonic  aspect  ratio  parameter 

20 

BA/K,  subsonic  aspect  ratio  parameter,  or  Bm,  supersonic 
sweep  parameter 
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TABLE  12.  PARAMETER  DATA  REGION  BS  FOR  SUBROUTINES 
USPANF  AND  BNLDSF  (CONa) 


Location 

Description 

BNLDSF 

BS(1)  is  equivalent  to  TC0M(2973) 

1 

Sg,  area  of  theoretical  surface  (ft^) 

2 

Xle*  fuselage  station  of  root  chord  leading  edge  (in.) 

3 

bg,  span  of  surface  (ft) 

4 

PZN/PZW1 , ratio  of  fuselage  nose  normal  load  to  wing 
normal  load  due  to  angle  of  attack 

5 

A,  balance  equation  factor 

6 

B,  balance  equation  factor 

7 

Cf  balance  equation  factor 

8 

Df  balance  equation  factor 

9 

E,  balance  equation  factor 

10 

F,  balance  equation  factor 

11 

W0P,  exposed  wing  panel  weight  per  side  (lb) 

12 

K(j,  flap  lift  effectiveness  parameter 

13 

PZW1,  wing  flexible  normal  airload  due  to  angle  of 
attack  (lb) 

14 

PZW2,  wing  rigid  normal  airload  due  to  flap  deflection  (lb) 

13 

PZW3,  wing  increment  flexible  normal  airload  due  to  flap 
deflection  (lb) 

16 

PZW4,  wing  flexible  normal  airload  due  to  vertical 
acceleration  (lb) 

17 

PZH,  horizontal  tailload  (lb) 

18 

airplane  mass  ratio 

19 

Kg,  gust  alleviation  factor 

20 

APZHQ,  incremental  horizontal  tailload  due  to  pitching 
acceleration  (lb) 
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TABLE  13.  THE  BD  DATA  REGION 


Location 

Description 

BD(1)  is  equivalent  to  TC0M(2993) 

1-4 

Array  (YA)  - subsonic  (BCj^/K)  table  values 

5-20 

Array  (YB)  - rigid  loading  table  values  (first) 

5-20 

Array  (YB)  - loading  at  weight  analysis  stations  (second) 

21-25 

Array  (YC)  - rigid  loading  values  at  aero  data  stations 

26-30 

Array  (ED)  - aero  data  stations  fraction  of  span 

31-43 

Array  (ESS)  - strip  boundary  stations  (first)  fraction 
of  span 

31-43 

Array  (ESS)  - weight  analysis  stations  (second) 
fraction  of  span 

44 

No.  of  wing  outer  panel  weight  strips  per  side 

45-57 

Array  (YSS)  - strip  running  loadings  (first) 

45-57 

Array  (YSS)  - torsional  loading  at  weight  analysis 
station  (second) 

58-68 

Array  (DPA)  - strip  alpha  loads 

69 

Flap  lift  center  of  pressure,  fraction  of  wing  chord  (X/C) 

70-80 

Array  (DPF)  - flap  effects  strip  loads 

81-88 

Subroutine  FC0DM2  factors 

81-84 

Array  (T) 

85-88 

Array  (YX) 

89 

PW0  - summated  strip  alpha  loads  to  centerline  (first) 

89 

Summated  AFW  (1,1)  - exposed  rigid  alpha  loads  (second) 
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TABLE  13.  THE  BD  DATA  REGION  (OONT) 


Location 

Description 

90 

*W1  " *W0  *ess  carry-over  strip  load  (first) 

90 

Simulated  AFWF  (1,1)  - exposed  flexible  alpha  loads 
(second) 

91 

Simulated  wing  outer  panel  weights  (first)  (lb) 

91 

Ratio  of  exposed  alpha  load  to  total  (second)  (dec)  a 1.0 

92 

Atmospheric  density  (Rho)  (slugs/feet  cubed) 

93 

Factor  (NB-I)  - number  of  strips  counter 

94 

Speed  of  Sound  (ah)  (knots) 

95-105 

Array  (ECP)  - strip  center  of  pressure,  fraction  of  span 

106-116 

Array  (ULN)  - running  load  at  weight  analysis  stations 

117-129 

Array  (DXS)  - swept  torsion  arms  (in.) 

130-140 

Array  (UTN)  - torsional  loading  at  weight  analysis 
station  (in.) 

141 

Parameter  K sweep 

142 

Fraction  of  swept  chord  of  reference  axis,  X/C 

143 

Fraction  of  swept  chord  of  center  of  pressure,  X/C 

144 

Swept  chord  (in.) 

145 

Expression:  1.0-Ks  (1.0- (X/C)  of  RA-(X/C)s  of  RA) 

146 

Integral  of  Array  (ULN)  versus  (Y/b/2)  - tip  to  root 

147 

KUZ,  span-loading  normalizing  parameter 

148 

Swept  bending  moment  at  body  side  (in. /lb) 

149 

Swept  torsion  moment  at  body  side  (in. /lb) 
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TABLE  13.  THE  BD  DATA  REGION  (CONCL) 


1 \ 


Location 

Description 

150 

Empennage  swept  torsion  arm  at  body  interface  (in.) 

151 

Body  nose  load  due  to  vertical  gust  (lb) 

152 

Exposed  wingload  due  to  vertical  gust  (lb) 

153 

Body  carry-over  load  due  to  vertical  gust  (lb) 

154 

Horizontal  tailload  due  to  vertical  gust  (lb) 

155 

Parameter  KPA 

156 

Wing  semispan  (first)  (in.) 

156 

Horizontal  tail  semispan  (second)  (in.) 

156 

Verticil  tail  span  (third)  (in.) 

157 

Wing  reference  axis  sweep  angle  (first)  (deg) 

157 

Horizontal  tail  reference  axis  sweep  angle  (second)  (deg) 

157 

Vertical  tail  reference  axis  sweep  angle  (third)  (deg) 

158 

Swept  distance  from  vertical  tail  root  to  horizontal 
tail  plane  (in.) 

159 

Gust  load  factor 

160 

Wing  inertial  aeroelastic  load  due  to  vertical  gust  (lb) 
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TABLE  14.  THE  BU  DATA  REGION 


Location 

Description 

BU(1)  is  equivalent  to  TcOM  (3153) 

1-60 

Data  for  wing  flexible  alpha  effects 

1 

Flexible  lift  curve  slope 

2 

X-distance  from  apex  to  total  load  CP  (in.) 

3 

Y-distance  from  CL  to  exposed  load  CP  (in.) 

4 

X-distance  from  apex  to  exposed  load  CP  (in.) 

5 

X-distance  from  apex  to  carry-over  load  CP  (in.) 

6 

Airload  unit  shear  at  side  of  body  (dec) 

7 

Airload  unit  bending  moment  at  side  of  body  (in.) 

8 

Airload  unit  torsion  moment  at  side  of  body  (in.) 

9-21 

13  load  stations,  percent  span  from  tip  to  £ (dec) 

22-34 

Airload  unit  siiears  from  tip  to  centerline  (dec) 

35-47 

Airload  unit  bending  moments  - swept  (in.) 

48-60 

Airload  unit  torsion  moments  - swept  (in.) 

61-107 

Data  for  wing  rigid  flap  effects 

61 

Parameter  KBF 

62 

X-distance  from  apex  to  total  load  CP  (in.) 

63 

Y-distance  from  CL  to  exposed  load  CP  (in.) 

64 

X-distance  from  apex  to  exposed  load  CP  (in.) 

65 

X-distance  from  apex  to  carry-over  load  CP  (in.) 

66 

Airload  unit  shear  at  side  of  body  (dec) 

67 

Airload  unit  bending  moment  at  side  of  body  (in.) 
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TABLE  14.  THE  BU  DATA  REGION  (OONT) 


Location 

Description 

68 

Airload  unit  torsion  moment  at  side  of  body  (in.) 

69-81 

Airload  unit  shears  from  tip  to  centerline  (dec) 

82-94 

Airload  unit  bending  moments  - swept  (in.) 

95-107 

Airload  unit  torsion  moments  - swept  (in.) 

108-166 

Data  for  horizontal  tail  effects 

108 

Horizontal  tail  lift  curve  slope 

109 

Parameter  KH(B) 

110 

YH  distance  from  centerline  to  horizontal  tailload  CP  (in.) 

111 

XH  distance  from  HT  apex  to  HT  load  CP  (in.) 

112 

Airload  unit  shear  at  side  of  body  (dec) 

113 

Airload  unit  bending  moment  at  side  of  body  (in.) 

114 

Airload  unit  torsion  moment  at  side  of  body  (in.) 

115-127 

Thirteen  HT  load  stations  - percent  spantip  to 
centerline  (dec) 

128-140 

Airload  unit  shears  from  tip  to  centerline  (dec) 

141-153 

Airload  unit  bending  moments  - swept  (in.) 

154-166 

Airload  unit  torsion  moments  - swept  (in.) 

167-225 

Data  for  vertical  tail  effects 

167 

Vertical  tail  lift  curve  slope 

168 

Parameter  KV(B) 

169 

Z-distance  from  VT  root  to  vertical  tailload  CP  (in.) 
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TABLE  14.  THE  BU  DATA  REGION  (QONT) 
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TABLE  14.  THE  BU  DATA  REGION  (CONCL) 


Location  Description 

365-377  Airload  unit  shears  from  tip  to  oenterline  (dec) 

378-390  Airload  unit  bending  moments  - swept  (in.) 

391-403  Airload  unit  torsion  moments  - swept  (in.) 

404-449  Data  for  inertial  aeroelastic  effects 

404  X-distance  from  apex  to  total  load  CP  (in.) 

405  Y-distance  from  centerline  to  exposed  load  CP  (in.) 

406  X-distance  from  apex  to  exposed  load  CP  (in.) 

407  X-distance  from  apex  to  carry-over  load  CP  (in.) 

408  Airload  unit  shear  at  side  of  body  (dec) 

409  Airload  unit  bending  at  side  of  body  (in.) 

410  Airload  unit  torsion  at  side  of  body  (in.) 

411-423  Airload  unit  shears  from  tip  to  centerline  (dec) 
424-436  Airload  unit  bending  moments  - swept  (in.) 

437-449  Airload  unit  torsion  moments  - swept  (in.) 

450-454  Data  for  subroutine  SPABMF 

450  Wing  total  flexible  alpha  airload  (PZW1)  (lb) 

451  Wing  total  rigid  flap  airload  (PZW2)  (lb) 

452  Wing  total  aeroelastic  flap  airload  (PZW3)  (lb) 

453  Wing  total  inertia  aeroelastic  airload  (PZW4)  (lb) 


455-494 

495-500 


Previously  inputed  El  and  GJ  arrays 
Are  not  used 


TABLE  15.  OUTPUT  B0  DATA  REGION 


Location 

Description 

B0(1)  through  B0(18O)  is  equivalent  to  TG0M(3653)  through  TC0M(3832) 

1 

Condition  number 

2 

Fuselage  ID  number 

3 

PZN>  body  nose  normal  load  (lb) 

4 

Pyjj,  body  nose  side  load  (lb) 

5 

Xj^,  body  nose  CP  (fus  sta)  (in.) 

6 

PzW(B)/2,  exposed  wing  panel  load  (per  side)  (lb) 

7 

%(B)»  spanwise  CP  of  exposed  wing  panel  load  (BL)  (in.) 

8 

X^B)»  exposed  wing  panel  CP  (fus  sta)  (in.) 

9 

PZB(W)’  k 0 ^ “rry-over  load  (lb) 

10 

XB(W)’  ^ carry-over  CP  (fus  sta)  (in.) 

11 

PZHr/2,  horizontal  tail  panel  load  (per  side)  (lb) 

12 

Y^,  spanwise  CP  of  horizontal  tail  panel  load  (BL)  (in.) 

13 

X^p  horizontal  tailload  CP  (fus  sta)  (in.) 

14 

Mxv(H)»  incremental  unsynmetrical  horizontal  tail  rolling 
moment  (for  T-tail  and  fus)  (in. /lb) 

15 

Pyvr»  vertical  tail  side  load 

16 

lyj,  vertical  tail  spanwise  CP  (WL)  (in.) 

17 

Xyp,  vertical  tail  CP  (fus  sta)  (in.) 

18 

N^,  airplane  normal  load  factor 

19 

Ny,  airplane  side  load  factor 

20 

• 2 
Q,  airplane  pitching  acceleration  (rad/ sec  ) 
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TABLE  15.  OUTPUT  B0  DATA  REGION  (GONT) 


Location 

Description 

21 

ft,  airplane  yawing  acceleration  (rad/sec2) 

22 

Wing  ID  number 

23 

Yjjiv*  wing-body  interface  station  (BL)  (in.) 

24 

SzBW»  wing  shear  load  at  side  of  body  station  (lb) 

25 

Mggty,  exposed  wing  panel  rolling  moment  at  side  of 
body  station  (in. /lb) 

26 

MygW,  exposed  wing  panel  pitching  moment  at  intersection 
of  load  reference  line  and  side  of  body  station  (in. /lb) 

Wing  shears  and  moments  at  stations  Y^A  (1)  through  YWA(12)  along  lo*d 
reference  line: 

27 

yWA  (1)*  first  wing  station  (in.) 

28 

Szw(l),  shear  at  station  Yj.,a(1)  (lb) 

29 

MxA(1),  bending  moment  at  station  Y^U)  (in. /lb) 

30 

MyA(D,  torsional  moment  at  station  Y^A(1)  (in. /lb) 

31  - 74 

Station,  shear,  and  moments  in  same  order  as  27  through 
30,  for  next  11  wing  stations 

75 

Horizontal  tail  ID  number 

76 

Yg^,  horizontal  tail-body  interface  station  (BL)  (in.) 

77 

SZBH*  horizontal  shear  at  side  of  body  station  (lb) 

78 

Mxbh*  exposed  horizontal  tail  panel  rolling  moment  at 
side  of  body  station  (in. /lb) 

79 

%BH»  exposed  horizontal  tail  panel  pitching  moment  at 
intersection  of  load  reference  line  and  side  of  body 
station  (in. /lb) 
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TABLE  15.  OUTPUT  B0  DATA  REGION  (GQNCL) 


Location 

Description 

Horizontal  tail  shears  and  moments  at  stations  Yha(1)  through  Y^(12) 
along  load  reference  line:. 

80 

Y„.(l),  first  horizontal  tail  station  (in.) 

HA 

81 

SZH,  shear  at  station  YHa(1)  (lb) 

82 

M^,  bending  moment  at  station  YHa(1)  (in. /lb) 

83 

MyA,  torsional  moment  at  station  YHa(1)  (in. /lb) 

84  - 127 

Station,  shear,  and  moments  in  same  order  as  80  through 
81,  for  next  11  horizontal  tail  stations 

128 

Vertical  tail  ID  number 

129 

Zgy,  vertical  distance  from  vertical  tail  root  chord 
station  to  vertical  tail-body  interface  station  (in.) 

130 

SYBV»  vertical  tail  shear  at  vertical  tail-body 
interface  station  (lb) 

131 

Mxbv»  exposed  vertical  tail  panel  rolling  moment  at 
vertical  tail-body  interface  station  (in. /lb) 

132 

MzbV»  exposed  vertical  tail  panel  yawing  moment  at  the 
intersection  of  the  load  reference  line  and  the  vertical 
tail -body  interface  station  (in. /lb) 

Vertical  tail  shears  and  moments  at  stations  Zva(1)  through  Zva(12) 
along  load  reference  line 

133 

ZVA(1),  first  vertical  tail  station  (in.) 

134 

Sz  (1) , shear  at  station  Z v (1)  (lb) 

135 

M^a(1),  bending  moment  at  station  Zy^(l)  (in. /lb) 

136 

M^l),  torsional  moment  at  station  2^(1)  (in. /lb) 

137  - 180 

Station,  shear,  and  moments  in  same  order  as  133  through 
136,  for  next  11  vertical  tail  stations 
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TABLE  16.  INTEGER  DATA  REGION  (ND) 

Description 


ND(1)  is  equivalent  to  TG0M(4201) 

1-12  Integers 

13-46  Control  factors  (refer  to  Table  1) 

47-99  Are  not  used 

100  BFCNTL  - Case  number  counter  (NC) 

101  Integer  Variable  (I) 

102  Integer  Variable  (J) 

103  USPANF  - N4EIGJ  - 20 

103  BFOTL  - Punched  card  index  number  col  2-12  (IC) 

104  Punched  card  ID  number  col  73-80  (IDC) 

105  BFOJTL,  USPANF  - SIC  matrix  calculation  option  (IACALCS) 

106  BFtm,  BNLDSF  - Condition  type  factor  (NI) 

107  BFCNTL,  RERDAT,  USPANF  - RERDAT  entry  factor  (IR) 

108  BFCNTL,  USPANF  - (NF  » -1) 

109-131  BFQJTL  - Condition  selection  array  (No  * -1,  Yes  - 1) 

132  RERDAT,  USPANF  - Integer  variable  (I) 

133  Integer  variable  (J) 

134  USPANF  - Integer  variable  (K) 

135  Integer  variable  (L) 

136  Integer  variable  (M) 

137  USPANF,  RERDAT,  BNLDSF,  SPABMF  - Wing  position  factor  (IP) 

(Refer  to  note) 
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TABLE  16.  INTEGER  DATA  REGION  (ND)  (OONT) 


Location 


Description 


138 

139 

140 

141 

142 

143 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 


FC0EM2  - Interpolation  factor  (N3) 
Interpolation  factor  (N4) 

Integer  variable  (I) 

Integer  variable  (K) 

Integer  variable  (L) 

O0DIM2  - Interpolation  factor  (Nl) 

Integer  variable  (J) 

Integer  variable  (JJ) 

Integer  variable  (K) 

Integer  variable  (L) 

Integer  variable  (M) 

USPANF  - Repeat  at  ive  operation  counter  (NT) 
BNLDSF  - Repeatative  operation  crunter  (NT) 
Integer  variable  (I) 

Integer  variable  (J) 

Integer  variable  (K) 

SPABMF  - Distribution  selector  (ID) 


155 

Integer  variable  (I) 

156 

Integer  variable  (J) 

157 

Integer  variable  (K) 
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TABLE  16.  INTEGER  DATA  REGION  (ND)  (CONCL) 


Location 

Description 

158 

Integer  variable  (L) 

159 

Integer  variable  (M) 

160-200 

Not  used 

NOTE:  (-1) 

■ fixed,  (0)  ■ VS  wing  fwd,  (1)  - VS  wing  aft 
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TABLE  17.  LISTING  OF  CONDITION  PUNCHED  OUTPUT  COMPONENT  LIMIT  AIRLOADS  DATA  FOR  SWEEP 
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TABLE  18.  CONDITION  ARRAY  ID  (CN) 
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CASE  7 

NEGATIVE  LfAO  FACTOR  BALANCED  FLIGHT  MANEUVER  AT  ALT(l,2) 
CASE  4 = MHU*ll 
CASE  5 = mh( 1,2) 

CASF  f HNs.qQ 


TABLE  18.  CONDITION  ARRAY  ID  (CN)  (CONT) 
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TABLE  18.  CONDITION  ARRAY  ID  (CN)  (CONCL) 
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TABLE  20.  INPUT  BF  DATA  SET 
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TABLE  21.  INTERPOLATED  AERODYNAMIC  DATA 

U B SONTLC  LV*  t A 
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TABLE  21.  INTERPOLATED  AERODYNAMIC  DATA  (CONT) 
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TABLE  21.  INTERPOLATED  AERODYNAMIC  DATA  rOONCL) 
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TABLE  22.  WFLEX  LOADING  MATRICES 
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TABLE  22.  WFLEX  LOADING  MATRICES  (CCNT) 


‘>1 

9 

<*i 

o 

N 

9 

IN 

9 

IN 

9 

IN 

9 

T 

4* 

T 

X 

T 

44 

T 

44 

T 

44 

T 

41 

X 

X 

i*i 

X 

X 

X 

o 

a 

9 

X 

o 

*4 

0 

X 

iN 

X 

•4 

X 

a 

X 

X 

X 

rf> 

X 

X 

X 

• 

t 

• 

• 

• 

• 

.*i 

X 

•4 

a 

•4 

» 

1 

i 

*4 

IN 

N 

IN 

IN 

3 

o 

o 

X 

9 

o 

X 

1 

1 

1 

T 

t 

44 

Ui 

X 

44 

Jd 

X 

•3 

X 

a 

X 

X 

P4 

X 

IN 

a 

Cl 

p4 

O 

a 

X 

X 

r* 

X 

p4 

i*i 

X 

04 

X 

* 

X 

X 

a 

a 

• 

• 

• 

• 

• 

• 

t 

P4 

m 

X 

r4 

t 

1 

i 

<N 

IN 

IN 

-N 

IN 

0 

1 

? 

f 

/*N 

f 

T 

f 

44 

X 

44 

K 

44 

44 

44 

1*1 

o 

o 

S 

<3 

X 

X 

*4 

(A 

.N 

K 

o 

p4 

£ 

IX 

N* 

X 

*n 

N 

X 

f» 

X 

* 

p4 

0> 

X 

t 

• 

o 

• 

$ 

in 

• 

r* 

• 

4 

• 

u> 

1*1 

X 

1 

U. 

Q 

f 

*4 

1 

IN 

IN 

O 

4 

9 

cx 

o 

■N 

9 

•N 

9 

9 

a 

T 

44 

T 

X 

T 

144 

4 

T 

44 

T 

44 

T 

44 

*4 

o 

. j 

X 

X 

«4 

♦ 

> 

x 

X 

X 

a 

■i. 

a 

o 

X 

X 

*4 

u> 

o 

o 

a 

!»• 

X 

1*1 

9 

9 

a 

SI 

X 

X 

X 

St 

• 

• 

t 

X 

• 

• 

• 

X 

x 

X 

f» 

P4 

111 

p4 

IN 

P4 

1 

aC 

1 

i 

ac 

t- 

rx 

.'U 

o4 

St 

IN 

IN 

•N 

St 

3 

9 

■3 

o 

O 

9 

T 

T 

1 

1 

1 

T 

2 

a* 

44 

X 

d 

44 

3 

44 

a 

<3 

9 

P* 

X 

p4 

3 

P* 

p4 

9 

p* 

X 

•4 

a 

r» 

*4 

a 

< 

irt 

X 

i*i 

o 

X 

a 

o 

X 

N 

X 

0 

i«i 

a 

-J 

• 

• 

t 

t 

• 

• 

r» 

r» 

*4 

X 

N 

x4 

#4 

«*4 

i 

44 

1 

i 

-i 

X 

X 

in 

p4 

p4 

u. 

N 

IN 

IN 

p4 

9 

9 

9 

O 

9 

9 

X 

T 

T 

T 

a 

1 

T 

T 

44 

UJ 

44 

44 

X 

44 

X 

-4 

x 

X 

X 

9 

p4 

*4 

a. 

9 

«4 

U4 

t 

p4 

<3 

® 

X 

a 

a 

X 

.<1 

a 

* 

O 

o 

9 

a> 

p4 

1*1 

<»1 

IN 

a 

o 

X 

• 

• 

• 

• 

• 

• 

X 

p4 

#4 

•4 

p4 

*4 

1 

i 

l 

<N  IN 

-4  4 

X I*) 

■N  IN 

«n  i*i 

l«!  1*1 

0 o 

1 1 

ff 

if 

ff 

?? 

f? 

uj  UJ 

44  44 

44  04 

44  44 

X X 

X X 

a x 

X * 

a x 

•N  1*1 

p4  p4 

X 

*4 

X X 

iN 

i»  a 

* 

o a 

*4 

9 x 

IN 

^ X 

1*1 

<*1  O 

IN  X 

X x 

X X 

x in 

X a 

X 1*1 

9 IN 

uS  a 

i*>  a 

a x 

X o 

X <3 

X x 

•A  o 

iN  1*1 

x a 

Xf» 

H 

• • 
CO  m 

N 

• • 

N 

• • 
r*  f 

ri 

• • 
* ^ 

N 

• • 

H 

t • 

• 

• 

• 

i t 

• 

1 • 

• 

• T 

• 

Z 

N <N 

8 

X -n 

8 

X «*1 

8 

1*1  iN 

8 

1*1  IN 

8 

1*1  in 

2 

?? 

if 

2 

ff 

2 

?f 

ff 

2 

ff 

2 

UJ  X 

2 

X X 

x 

X X 

i 

X X 

V 

X X 

Z 

X X 

3 

O IN 

3 

X X 

5 

X O 

a 

X*  IN 

1 3 

O 9 

a 

^ o 

■j 

X > 

r»  X 

X 

X 9 

X 

X*  X 

X 

1*1  X 

X 

1*1  X 

O 

X X 

s 

9 9 

o 

^ *4 

o 

N X 

a 

2 O' 

o 

a x 

‘J 

5 9 

u 

■4  p4 

u 

• r» 

u 

x a 

o 

O 9 

o 

• 9 

N X 

x > 

IN  (A 

9 N 

f»  (*1 

» n ^ 

• • 

• • 

• • 

• • 

• • 

• • 

9 IN 

X IN 

«n  m 

X X 

1*1  x 

X «N 

i 1 

1 1 

1 1 

128 


TABLE  22.  WFLEX  LOADING  MATRICES  (OONCL) 
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TABLE  23.  SPANWISE  LOADING  DISTRIBUTION 
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TABLE  23.  SPANWISE  LOADING  DISTRIBUTION  (CONT) 
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TABLE  23.  S PANWISE  LOADING  DISTRIBUTION  (GONCL) 
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TABLE  24.  USPANF  SPANWISE  UNIT  DISTRIBUTIONS 
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TABLE  24.  USPANF  SPANWISE  UNIT  DISTRIBUTIONS  (CONT) 
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TABLE  24.  USPANF  S PANWISE  UNIT  DISTRIBUTIONS  (CONT) 
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TABLE  24.  USPANF  SPANWISE  UNIT  DISTRIBUTIONS  (CONT) 
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TABLE  24.  USPANF  SPANWISE  UNIT  DISTRIBUTIONS  (CCNCL) 
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TABLE  25.  COMPONENT  LIMIT  AIRLOADS  AND  CENTERS  OF  PRESSURE 
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TABLE  26.  WING  LIMIT  AIRLOAD  DISTRIBUTIONS 
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TABLE  27.  HORIZONTAL  TAIL  LIMIT  AIRLOAD  DISTRIBUTIONS 
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TABLE  33.  STRUCTURAL  INFLUENCE  COEFFICIENTS  (IN. /LB)  MATRIX  (CONCL) 
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TABLE  34.  DTHEDA  MATRIX 
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TABLE  34.  DTI  EDA  MATRIX  (CONCL) 
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TABLE  35.  SICRAR  MATRIX  (RAD/LB) 
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TABLE  35.  SICRAR  MATRIX  (RAD/LB)  (CONCL) 
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TABLE  36.  DATA  FROM  SUBROUTINE  WFLEX 
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Section  IV 


PROGRAM  USAGE 


GENERAL 


The  flexible  airloads  stand-alone  program,  BFCNTL,  is  designed  to 
perform  functions  similar  to  the  airloads  module,  BLCNTL,  of  SWEEP;  that  is, 
the  program  computes  airloads  on  the  airplane  components  for  selected  flight 
conditions.  Punched  card  input  variable  data  are  generated  by  SWEEP.  Opera- 
tion of  the  program  requires  the  user  to  prepare  a program  control  factors 
card  and  to  set  up  the  program  decks.  Punched  card  output  data  are  in  a 
format  that  is  compatible  for  use  as  an  optional  input  to  the  SWEEP  program. 


PROGRAM  DECK  SETUP 

The  program  deck  setup  is  illustrai  ed  in  Figures  6 and  7.  The  total 
deck  setup  must  follow  the  blocked  order  shown.  Twelve  subordinate  subrou- 
tines may  be  arranged  in  any  order,  but  the  total  subroutine  block  must  be 
immediately  behind  the  main  program  deck.  When  only  one  case  is  to  be  run, 
the  execute  card  must  follow  the  last  card  of  the  DP  data  set. 

When  multiple  cases  are  to  be  run,  subsequent -case  data  must  be 
arranged  as  shown  in  Figure  7 and  placed  immediately  behind  the  first  case. 
The  execute  card  must  then  follow  the  last-case  DP  data  termination  card. 


DATA  INPUT  DESCRIPTION 


An  ND  control  factors  card  must  be  provided  for  each  case  data  block  of 
a multiple-case  computer  run.  It  must  be  the  first  card  of  each  case  data 
block.  The  format  of  this  card  is  discussed  on  page  , and  illustrated 
in  Table  1. 


Data  set  BC  follows  the  ND  card.  Data  set  BF  follows  data  set  BC.  After 
data  sets  BC  and  BF  are  read  into  core , subroutine  RERDAT  merges  the  data  and 
foims  a new  BC  data  set.  The  initial  BF  data  set  is  not  disturbed.  It  is 
necessary,  therefore,  to  repeat  the  entire  input  BC  data  set  for  each 
subsequent -case  data  block  of  a multiple-case  computer  run.  Data  sets  DT,  DB, 
DF,  and  DP  follow  data  set  BF  in  the  order  indicated  in  Figure  7. 
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All  of  the  variables  in  data  sets  BC  through  DP  are  in  the  form  to  be 
read  by  subroutine  DECRD.  This  subroutine  and  data  card  format  are  described 
in  Section  III.  Data  set  BF  for  subsequent -case  data  blocks  of  a multiple- 
case  computer  run  needs  to  contain  only  appropriate  change  variables  with  a 
DECRD  termination  punch  in  Column  1 of  the  last  card  in  the  data  set,  or  a 
DECRD  termination  card  if  there  are  to  be  no  changes.  Data  sets  DT,  DB,  DF, 
and  DP  for  subsequent -case  data  blocks  require  only  the  DECRD  termination 
card  for  each  set,  since  these  data  sets  are  read  into  core  for  the  first 
case  and  remain  undisturbed. 


PROGRAM  CONTROL  FACTORS  AND  OPTIONS 

The  control  of  the  program  to  obtain  the  user's  desired  options  is 
accomplished  by  the  input  data  contained  on  the  ND  control  factors  card.  This 
card  is  prepared  by  the  program  user  and  must  be  prepared  with  a great  deal  of 
care  to  obtain  desired  output  data. 

A complete  description  of  the  data,  locations  by  card  column,  and  the 
options  available  are  shown  in  Table  1.  All  of  the  variables  in  the  list  are 
in  integer  form  and  are  entered  in  two-column  fields.  Single-character 
variables  are  entered  in  the  second  column  of  the  field. 

Factor  ND(13)  indentifies  the  air  vehicle  class.  Factors  ND(14)  and 
ND(15)  identify  the  wing  type  (fixed  or  variable  sweep)  and  the  vertical  tail 
type  (single,  dual,  or  T-tail). 

Factors  ND(23)  through  ND(27)  instruct  the  program  to  compute  loads  on 
selected  airframe  components  or  on  all  airframe  components. 

Factors  ND(28)  through  ND(36)  provide  for  the  selection  of  the  types  of 
conditions  to  be  computed.  Figure  1 and  Table  2 show  the  types  of  conditions 
that  may  be  selected  for  each  of  the  11  case  numbers  available  for  entry  in 
ND(46) . Any  condition  type  not  specifically  defined  by  the  chosen  case 
number,  ND(46),  will  not  be  computed,  even  if  the  type  factors  ND(28)  through 
ND(36)  are  entered  as  "yes." 

Factor  ND(40)  is  used  to  select  or  reject  the  punched  card  output.  Since 
the  primary  purpose  of  the  program  is  to  obtain  the  punched  output  for  use  as 
an  optional  input  to  the  SWEEP  program,  it  would  be  expected  that  the  user 
would  always  put  a "yes"  in  ND(40) . The  factor  ND(41)  provides  for  a print- 
out listing  of  the  data  contained  in  the  punched  card  output. 

Factors  ND(44)  and  ND(45)  provide  for  the  optional  printout  of  FLXSIC 
and  WFLEX  matrices. 
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Factor  ND(42)  defines  the  number  of  stations  and  El  and  GJ  values  that 
are  used  by  subroutines  FLXSIC  and  WFLEX  to  describe  wing  stiffness  distribu- 
tion. The  ND(42)  entry  should  always  be  20. 

Factor  ND(43)  defines  the  number  of  chordwise  loads  strips  along  the 
exposed  wing  which  are  used  in  the  calculation  of  the  aeroelastic  loadings. 
The  ND(43)  entry  should  always  be  10. 

Card  columns  7 through  20  and  49  through  54  should  be  left  blank. 

BC  DATA  SET  USAGE 

A complete  description  of  input  BC  data  set  is  presented  in  Table  3. 

All  of  the  variables  in  the  list  of  inputs  are  in  the  form  to  be  read  by 
subroutine  DECRD.  Expanded  explanations  for  those  variables  which  might  be 
subject  to  misinterpretation  are  in  the  following  paragraphs. 

The  following  data  locations  need  be  entered  only  if  a variable  sweep 
wing  in  the  forward  position  is  to  be  analyzed.  Applicable  case  numbers  for 
this  configuration  are  8,  9,  10,  11,  as  indicated  in  Table  2. 

BC(2) 

BC(5) 

BC(7) 

BC(9) 

BC(25)  through  BC(30) 

BC(69)  through  BC(88) 

If  a fixed  wing  or  a variable  sweep  wing  in  the  aft  position  is  to  be 
analyzed,  these  input  data  locations  may  be  eliminated  or,  preferably  entered 
as  zeros. 

The  following  input  data  locations  are  not  used  for  any  input  data  and 
may  be  eliminated  or,  preferably  entered  as  zeros. 

BC(58)  through  BC(68) 

BC(89)  through  BC(99) 

BC(126)  through  BC(136) 

BC(158)  through  BC(165) 

A separate  BC  data  set  must  be  entered  for  each  subsequent-case  data 
block  of  a multiple-case  computer  run. 
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BF  DATA  SET  USAGE 


A complete  description  of  input  BF  data  set  is  presented  in  Table  4.  All 
of  the  variables  in  the  list  of  inputs  are  in  the  fora  to  be  read  by  subrou- 
tine DECRD.  Expanded  explanations  for  those  variables  which  might  be  subject 
to  misinterpretation  are  in  the  following  paragraphs. 

BF(2)  and  BF(25)  through  BF(46)  data,  locations  need  be  entered  only  if  a 
variable  sweep  wing  in  the  forward  position  is  to  be  analyzed.  The  applicable 
case  numbers  for  this  configuration  are  8,  9,  10,  and  11,  as  indicated  in 
Table  2. 

If  a fixed  wing  or  a variable  sweep  wing  in  the  aft  position  is  to  be 
analyzed,  these  input  data  locations  are  eliminated  or,  preferably  entered 
as  zeros. 

Data  locations  BF(5)  through  BF(14)  and/or  BF(27)  through  BF(36)  contain 
the  dead  weights  of  the  10  chordwise  strips  along  the  exposed  wing  span, 
starting  with  the  most  inboard  strip  and  ending  with  the  most  outboard  strip 
as  shown  in  Figure  2.  Data  locations  BF(15)  through  BF(24)  and/or  BF(37) 
through  BF(46)  contain  the  corresponding  chordwise  (X/C)  center-of-gravity 
values . 

Stations  along  the  elastic  axis  entered  in  BF(47)  through  BF(66),  and 
the  El  and  GJ  data  entered  in  BF(67)  through  BF(86)  and  BF(87)  through  BF(106) , 
should  contain  values  for  the  body  side  station,  a station  slightly  outboard 
of  the  body  side,  and  a station  at  or  near  the  tip.  Inclusion  of  these 
stations  insures  an  accurate  curve  fit  operation  by  subroutine  C0DIM2.  A 
maximum  of  20  stations  can  be  used. 


SEMIPERMANENT  DATA  SETS 

DT,  DB,  DF,  and  DP  data  sets  contain  all  of  the  required  aerodynamic 
data  required  to  develop  the  lifting  surface  load  distributions.  These  data 
are  semipermanent  (fixed  data)  and  are  the  same  as  the  data  contained  in 
the  SWEEP  program  data  bank.  These  data  are  described  in  Section  III  and  in 
Tables  5 through  8.  Data  cards  are  in  the  format  to  be  read  by  subroutine 
DECRD. 


OUTPUT  DATA  DESCRIPTION 

Primary  output  of  the  program  is  a deck  of  cards  containing  airplane 
component  limit  airloads,  centers  of  pressure  locations,  and  wing  and 
empennage  limit  airload  shears  and  moments  at  spanwise  stations  along  the 
load  reference  line  selected  for  weight  analysis.  The  punched  card  output  is 
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a format  for  direct  use  as  an  on  optional  external  ail  \oads  input  to  SWEEP. 
The  printout  of  other  data  is  also  included  for  visual  inspection  of  the 
results  of  final  and  intermediate  calculations. 


PUNCHED  CARD  OUTPUT 

Output  punched  cards  list  the  limit  airload  data  in  an  E- format  to  be 
read  by  subroutine  DECRD  in  SWEEP.  A sample  printout  of  the  data  contained 
on  the  punched  cards  is  shown  in  Table  17.  Each  card  contains  a DECRD  array 
index  nunber,  followed  by  five  decimal  data  items  and  the  card  identification 
nunber.  The  output  data  array  is  described  in  Table  15.  The  index  number 
punched  on  each  card  corresponds  to  the  location  within  the  array  for  the 
first  item  of  decimal  data  on  the  card.  A "yes"  in  ND  control  factor  ND(40) 
will  provide  the  punched  output  data  deck. 

Condition  identification,  the  component  loads  identification,  and  the 
punched  card  identification  numbers  are  assigned  by  the  program  and  are 
included  in  the  punched  output  deck  for  each  condition.  Component  loads 
identification  numbers  also  appear  in  the  printed  output.  These  identifica- 
tion numbers  described  and  located  under  "Punched  Card  Output,"  in  Section  III. 


PRINTED  OUTPUT  DATA 

Refer  to  Section  III  for  a detail  description  of  the  print  options  and 
the  printed  data  output,  including  the  data  of  primary  interest,  the 
intermediate -step  and  diagnostic  data,  and  sample  printouts. 
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Appendix  A 

PROGRAM  FLOW  CHARTS  AND  FORTRAN  LISTINGS 


FLOWCHART  USflGE 

The  automatically  generated  computer  program  flow  charts  (AUTOFLOW) 
presented  in  this  document  include  a table  of  contents,  flow  charts,  and 
FORTRAN  listings  of  all  routines  in  the  module.  Ihe  80-column  card  listings 
are  sequenced  and  grouped  by  routine. 


CROSS-REFERENCE  LIST 

Ihe  AUTOFLOW  table  of  contents  which  precedes  the  flow  charts  and  FORTRAN 
lists  serves  to  cross  reference  the  latter  two.  This  table  lists  the  following 
from  left  to  right: 

* The  card  identification  from  columns  73  through  80  of  this  card,  or  card 
sequence  number.  When  sequence  number  is  used  in  place  of  card  identi- 
fication, it  is  enclosed  in  parentheses. 

* Ihe  page  and  box  number  where  this  card  is  displayed  in  a flow  chart. 

* Ihe  FORTRAN  statement  number  from  colunns  1 through  S of  this  card. 

* Ihe  card  identification  (s)  or  sequence  number  (s)  of  the  card(s)  refer- 
ring to  this  card  (repeated  as  required). 

* Ihe  pages  and  box  numbers  where  the  cards  referring  to  this  card  are 
displayed  in  a flow  chart  (repeated  as  required). 


FLOW  CHARTS 

Ihe  flow  charts  produced  by  AUTOFLOW  use  USASI  conventional  symbols. 
Since  the  flow  charts  are  mechanically  drawn  from  the  program  source  deck, 
there  are  no  omissions  or  vague  generalizations  about  the  processing  within 
the  boxes. 
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Every  box  on  each  page  is  uniquely  numbered  and  may  be  referred  to  from 
elsewhere  in  the  program.  The  source  of  a reference  to  a box  will  be  indicated 
by  showing  the  page  and  box  number.  If  the  number  is  followed  by  an  asterisk, 
there  are  multiple  references  to  this  point,  and  the  others  may  be  found  by 
using  the  cross-reference  list. 


The  most-often-used  symbol  is  the  decision  box.  Like  all  boxes,  its  box 
number  is  above  and  to  the  right  of  the  box.  Its  FORTRAN  statement  nunber  is 
above  and  to  the  left  of  the  box.  The  decision  choices  for  the  paths  are 
printed. 


The  unconditional  transfer  connecter  has  its  page  number  destination 
printed  above  or  to  the  left  of  the  box  nunber  destination  within  the  connecter. 
If  there  is  a FORTRAN  statement  number  at  the  destination,  it  is  printed  below 
the  connecter. 


The  exit  box  example  shows  a connecter  from  page  9,  box  15. 


The  subroutine  call  box  includes  the  calling  sequence.  The  page  and  box 
numbers  of  the  flow  chart  of  the  called  subroutine  are  shown  on  the  left-hand 
side  of  the  box.  The  page  number  is  above  the  box  number. 


r 

i 

' i 

|| 
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The  note  box  encloses  comments  of  a functional  nature, 


as  differentiated  from  the  21  column  comments,  which  a *e  left  justified  without 
a box,  that  show  the  comment  cards  included  in  the  FORTRAN  deck. 


The  process  box  is  used  to  enclose  FORTRAN  arithmetic  statements. 


Input  and  output  are  shown  as  comnunicating  with  a device.  The  list  used 
follows,  if  appropriate: 


The  computed  G0  T0  becomes  a branch  table  showing  the  page  and  box  number 
of  each  of  the  ordered  branches. 


The  column  connecters  and  initial  connecters  are  the  only  boxes  without 
external  box  numbers.  The  function  of  the  initial  connecter  is  always  clear, 
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but  the  label  given  is  the  symbol  in  the  next  FORTRAN  card,  which  is  often 
blank. 


The  column  connecter  identifies  the  page  and  box  number  to  which  it 
connects . 


DIRECTORY  OF  FLOW  CHARTS  AND  FORTRAN  LISTINGS 


Page  Number  for: 

Program  or  Subroutine 

Flow  Charts 

Listings 

Table  of  Contents  and  References 

173 

BFCNTL 

183 

248 

RERDAT 

195 

251 

nncRD 

200 

253 

USPANF 

202 

253 

ATM0S 

219 

261 

FC0EM2 

221 

262 

C0DIM2 

224 

262 

WFLEX 

228 

264 

FLXSIC 

231 

265 

GLSQ 

235 
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